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PROJECT SUMMARY 


This volume, prepared by Technology Group, contains the final report on 
the satellite-based mobile communications project. The research effort was 
based on the comprehensive study of a class of continuous-phase modulated 
(CPM) signals used in conjunction with trellis-coded modulation (TCM), 
interleaving/interlacing, and coherent /differential/ noncoherent detection. 

The general framework of our tasks was to provide a detailed analysis of 
a bandwidth-efficient coding and modulation, and select a candidate system 
which will be capable of transmitting 4800 bits/s over a channel with a 5- 
KHz bandwidth subject to Rician fading. In particular, we were expected 
to provide: 

• A detailed analysis and design of a trellis-coded CPM transmission 
scheme with suboptimum detection, which could offer implementa- 
tional simplicity and the provision for interleaving/deinterleaving to 
combat the fading effects. 

• Development of software for the computation of power spectrum. 

• Design and development of a Doppler frequency shift compensator. 

• Design and development of a symbol timing recovery circuit. 

• Extensive simulations to be used for system parameter optimization. 

• Breadboard implementation of the transmission system. 

As outlined in our initial Quarterly Report, the project was organized 
into a set of tcisks and subtasks based on which, a comprehensive exami- 
nation of adl the issues was carried out. Accordingly, the following tasks 
were completed: extensive analysis and substantial simulations of partial 
response emd full response CPM signals, generalized precoding, trellis codes, 
interleaving and interlacing, suboptimum detection schemes, theoretical 
performance bounds, and filtersi development of a class of Doppler estima- 
tors; design of symbol synchronizer; extensive examinations the spectrum 
of the various CPM signals ; and, a comprehensive system optimization. 
Additionally, the hardware implementation of the final candidate system 


111 




has been undertaken^. In an effort to provide a detailed description of 
our work, a set of comprehensive and illustrative Quarterly Reports was 
produced and submitted. 

Our study has demonstrated that the synergy between TCM, which im- 
proves error probability, and CPM signals, which provide constant envelope 
and low spectral occupancy, provides a satisfactory solution to the problem 
of transmitting on mobile satellite channels. 

In principle, two implementations of this idea are feasible. The first one 
takes advantage of both the bandwidth efficiency and the power efficiency of 
CPM codes, by using a receiver wHch combines the trellis structure of TCM 
and that of CPM. In this situation, the complexity can grow very large and 
become qmckly unmanageable, so that a suboptimum solution should be 
devised. Moreover, the introduction of an interleaver/deinterleaver pair, 
which is known to be beneficicd fro transmission on fading channels, would 
not be allowed. To obviate the associated problems, we avoided maximum- 
likelihood decoding of the CPM signals, and instead, differentially demod- 
ulated the sign2ds symbol-by-symbol. By doing this, the power efficiency 
of CPM is not exploited, however, its spectral properties are retained, and 
interleaving/deinterleaving is made possible. 

This tr ansmi ssion scheme was extensively simulated under several con- 
ditions; full response and partial response formats, several frequency pulse 
shapes, with and without precoding, detection points, receiver filter shapes, 
TCM schemes, interleaving/deinterleaving sizes and depths, and various 
fading channel parameters were considered. 

The results of our studies show that, after a careful selection of the 
system parameters, on fading channels differentially-detected CPM offers 
an error performance which is essentially the same as differentially coherent 
PSK. Now, since PSK does not use the bandwidth in a very efficient way, 
it has to be band-pass filtered to meet the requirements of closely- spaced 
mobile-radio channels. As a result, its envelope is no longer constant, and 
its performance would be degraded by power amplifiers operated at or nezu 
saturation for better power efficiency. On the other hand, CPM with the 
parameters chosen is bandwidth efficient, so that it does not require narrow 
filtering and consequently offers constant or near-constant envelope. Hence, 

^Additional effort is required to finalize the hardware. 
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under the same mobile system configurations, the trellis coded CPM appears 
to offer a better performance than the filtered trellis coded 8DPSK. 
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1 Introduction and motivation of the work 

In this Report we consider a satellite-based mobile communication scheme 
bzised on continuous-phase modulated (CPM) signals used in conjimction 

with trellis-coded modulation (TCM). 

In satellite-based land mobile communication systems both bandwidth 
and power are Umited resources. In fact, these systems employ frequency- 
division multiple access (FDMA) with a given channel spacing (say, 5 KHz), 
and the fraction of out-of-band power should be very small to prevent in- 
terferences to adjacent channels. On the other hand, the satellite distance 
from earth, its power limitations, and the need for low-cost (and hence 
low-gain) mobile antennas puts a serious hmit on the power resources. Ad- 
ditionally, the fading environment of mobile communication further limits 
the power efficiency of the system. 

In a bandwidth- and power-limited environment, a bandwidth- and 
power-efficient coding/modulation scheme must be used. Trellis coded 
modulation (TCM) offers an attractive scheme. It combines the choice 
of a higher-order modulation scheme with that of a convolutional code, 
while the receiver, instead of performing demodulation and decoding in 
two separate steps, combines the two operations into one. As a result, the 
reliability of a digital transmission system is increased without increasing 
the transmitted power nor the required bandwidth. By using TCM, simple 
schemes can be designed that achieve significant power gains (from 3 to 6 
dB) without any bandwidth expansion. 

Due to the strictly bandlimited environment created by the mobile satel- 
lite channel, the signals to be used in conjunction with trellis codes must be 
chosen carefully. Besides having excellent spectral characteristics, the sig- 
nalg used should have constant envelope if nonlinear amplifiers are used for 
better power efficiency. A class of bandwidth-efficient signals that satisfies 
both constraints is offered by continuous-phase modulated (CPM) signals, 
based on phase modulation where phase continuity is introduced to reduce 
the bandwidth occupancy. 

The synergy between TCM, which improves error probability, and CPM 
signals, which provide constant envelope and low spectral occupancy, is 
expected to provide a satisfactory solution to the problem of transmitting 
on mobile satellite channels. 
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In principle, two implementations of this idea are feasible. The first one 
takes advantage of both the bandwidth efficiency and the power efficiency 
of CPM codes, by using a receiver which combines the trellis structure 
of TCM and that of CPM. In this situation, TCM and CPM can even 
be integrated in a single entity, but the number of states necessary for a 
trellis representation of these signals is given by the number of TCM states 
times the number of CPM states. As this number can grow very large, 
the complexity of the receiver may become quickly unmanageable, and a 
suboptimum solution should be devised. Moreover, the introduction of an 
interleaver/deinterleaver pair, which is well-known to introduce beneficial 
effects in digital transmission over channels affected by fading, would not 
be allowed without destroying the spectral properties of CPM. 

yVc can trade a decrease in complexity for a decrease in power efficiency 
(but not in bandwidth efficiency) by giving up maximum-likelihood decod- 
ing of the CPM signals, which are instead demodulated symbol-by-symbol. 
By doing this, we still take advantage of the spectral properties of CPM, 
and the introduction of an interleaving/ deinterleaver pair is now possible. 
This is an added attractive feature of this suboptimum solution, since it in- 
creases the protection of the transmitted signal from the effects of selective 
fading. 

We have considered three solutions: coherent, differentially-coherent, 
and non-coherent symbol-by-symbol detection. Since differentially coherent 
detection appeared to us to be the most promising technique, our work was 
based on a CPM/TCM system based on this kind of detection. This report 
describes our analysis, simulation, and implementation of differentially - 
coherent symbol-by-symbol detection of trellis-encoded continuous-phase mod- 
ulated signals. 

2 Channel models, CPM, and TCM 

Besides additive Gaussian noise, which is the standard environment for the 
analysis of coding schemes for the transmission of digital data or speech, 
there is a number of additional sources of performance degradation that 
must be taken into account to assess the merits of a proposed transmission 
scheme for mobile satellite channels. The most important among them are 
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• Doppler shifts. These are due to mobile vehicle motion. If differ- 
ential detection is used, the information-bearing phase turns out to 
be shifted by Jin amount 27r/(fT,, where l/T, is the data symbol rate 
and fd is the Doppler frequency shift, which for operation at L-band 
can be expected to be up to 200 Hz. At a rate of 2400 symbols per 
second the corresponding phase shift is 30°. 

• Fading and shadowing. The transmitted radio signal reaches the 
receiver through different paths caused by reflections from obstacles, 
yielding a signal whose components, having different phases and am- 
plitudes, may either reinforce or cancel each other. Shadowing is 
caused by the obstruction of radio waves by buildings, trees, and 
hills. 

• Adjacent channel interference. The 5-KHz mobile channel used 
for trjinsmission operates in a channelized environment. As a result, 
signals suffer from interference from signals occupying adjacent chan- 
nels. 

• Channel nonlinearities. Primarily because of the high-power am- 
plifier in the transmitter, operated at or near saturation for better 
power efficiency, the channel is inherently nonlinear. 

• Finite interleaving depth. In order to break up the error bursts 
caused by amplitudes fades of duration greater than symbol time, 
encoded symbols should be interleaved. Now, infinite interleaving 
provides a memoryless channel, but in practice the interleaving frame 
must be limited. In fact, for speech transmission the total coding/ de- 
coding delay must be kept below 60 ms in order not to cause percep- 
tual annoyance. If the depth of interleaving cannot be larger than 
the maximum fade duration anticipated, this causes a performance 
degradation. 
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2.1 Continuous Phase Modulation (CPM) 

A continuous-phase modulated signal is defined by 


s{t, a) = y-^cos(27r/ot + 0{t, a)) (1) 

where E, is the symbol energy, T, is the symbol time, and fo the carrier 
frequency. The trzmsmitted information is contained in the phase 

OO 

6{t,si) = 2nh 53 anq{t-nT,) (2) 

n=— OO 


with q{i) the phase-shaping pulse given by 

q{t) = / (3) 

^(t) is the frequency pulse with finite duration, i.e., g{t) is nonzero only 
for 0 < t < LTg, L the pulse length. The value of L contributes to the 
taxonomy of CPM signals, nzimely: 

• L = 1: Full-response signals 

• T > 1: Partial-response signals. 

It is common to assume 

[ = 1 / 2 - 


In (2), 

a = • • ■ , <X- 2 ) ^- 1 ? j ■ ■ ' 

denotes the symbol sequence sent to the CPM modulator. The symbols 
take values ±1, ±3, ■ • ■ , ±{M - 1), where M = 2”*, m a positive integer. 
The paraimeter h is called the modulation index. 
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2.2 Combining TCM With CPM 

If CPM signals are combined with an external convolutional encoder, or 
equivalently they are used as the signal constellation to be used in a trellis- 
coded modulation (TCM) scheme, a further improvement can be obtained. 
This new scheme is obtained by observing that at the output of the trellis 
encoder we get a multilevel signal which in turn can be used as the input 
to the continuous-phase modulator. The design of the coding scheme and 
of the modulator scheme should be performed jointly in order to maximize 
the Euclidean distance residting from the combination of the two. 

To implement TCM/CPM, we may want to take advantage of both the 
bandwidth efficiency and the power efficiency of CPM codes. However, 
the complexity of the resulting optimum demodulator can grow very large, 
so that a suboptimum solution should be sought, which trades a decrease 
in complexity for a decrease in power efficiency (but not in bandwidth 
efficiency). This is obtained by uncoupling the demodulation of CPM from 
decoding of TCM. 

To keep the demodulation problem separate from the decoding prob- 
lem, we consider estimating the CPM phases symbol-by-symbol, and using 
these estimates to build up a metric for the TCM decoder. Then, we use an 
optimum ( Viterbi) algorithm for decoding, the complexity of it being that 
of the TCM scheme only. This procedure obviously entails a loss of opti- 
mality, which is traded against a manageable receiver structure. However, 
as mentioned before, the spectral properties of CPM are preserved. 

We have first analyzed coherent and noncoherent symbol-by-symbol de- 
modulation of CPM signals. The result was that noncoherent demodulation 
entails a very large loss of optimality, while coherent demodulation, which 
offers a better performance, may still be too complex to implement. For 
these reasons, we decided to consider in more depth differential detection, 
which has a complexity comparable to noncoherent detection without its 
performance penalty, 

2.3 Differential detection of CPM 

The complex envelope of the received signal, say f(<) is delayed by T, 
seconds, transformed into its conjugate f*{t - T,), and multiplied by itself. 
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Then the real 8uid imaginary components of the signal f{t)f*{t — T,) are 
sampled every T, seconds. As a result, a discrete signal is obtained whose 

phase is 

2irfoTs + A^n + Vn 

where A^n represents the change over one symbol interval of the signal 
phase, cind rjn represents the chzmge in phase due to the noise. Under the 
assumption that /oT, is an integer number, estimate of this phase provides 
a noisy estimate of A9„, which is used to recover the information sequence. 

For full-response CPM (the case of precoded partial response was also 
considered) we have 

SO that the maximum value taken on by the phase shift is 


= irMM - 1). 

This quantity must be lower than ir because the phases are observed modulo 
2ir. Thus, we must choose 


2.4 Doppler phase shift removal 

Let us write the observed signal r(f) in the form 

r(t) = (5) 

where P(<) and u{t) are the amplitude fluctuation and the phase fluctua- 
tion, respectively. They accovmt for fading, noise, and intersymbol inter- 
ference. fd is the Doppler frequency shift. 

For differential processing, the signal (5) is synchronously sampled at 
times in = nT., n = 0, 1, . . . , and the following sequence is formed: 

( 6 ) 

where (f>d = 2-trfdT, is the Doppler phase shift, 

A^n = = TThOn 
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represents the phase fluctuations due to the transmitted data, and tj^ — 
i/(t„) - u{in-i)- FinaUy, 

Pn = 

(Perfect symbol synchronization is assumed.) 

It is seen from (6) that the presence of a Doppler frequency shift causes 
the phase of the sequence r„ to be altered by a term (f>d, added to the 
information sequence A^„. This has to be removed prior to demodulation. 

Two constraints are associated with the design of a Doppler compen- 
sation circmt for continuous-phase modtalated signals with application to 
mobile satellite communications; 

• Fast frequency acquisition is required. In fact, if data are transmit- 
ted in short bursts or packets, the acquisition time should not be a 
signiflcant portion of the burst interval. 

• Since multipath fading affects the propagation, the Doppler estimator 
must be insensitive to the signal amplitude fluctuations caused by 
fading. 

As a consequence of the first requirement, we have consider open-loop 
estimation structures. Three estimators, that trade robustness to Doppler 
shifts for complexity, were proposed and their performance analyzed. 


2.5 Symbol synchronization 

The operation and performance of a symbol synchronizer suitable for our 
transmission scheme have been analyzed and simulated. The idea imderly- 
ing the synchronizer scheme is the production of spectr^J lines in the CPM 
signal, as generated by passing it in a non-linear device. These spectral 
lines, or periodic-like contributions, provide a mechamsm for extracting 
symbol timing. In fact, the spacing of the lines is 1/T„ where T, is the 
symbol duration. 


3 Simulation results and conclusions 

The transmission scheme described here was extensively simulated under 
several conditions. Several frequency pulse shapes, receiver filter shapes, 
TCM schemes, interleaving/deinterleaving sizes and depths, and fading 
channels were considered. The results show that, for a careful selection 
of system parameters, on fading channels differentially-detected CPM of- 
fers an error performance which is essentially the same as differentially 
coherent PSK. Now, since PSK does not use bandwidth in a very efficient 
way, it has to be band-pass filtered to meet the requirements of closely- spaced 
mobile-radio channelization. As a result, its envelope is not constant, and 
its performance would be degraded by power amplifiers operated at or near 
saturation for better power efficiency. On the other hand, CPM with the 
parameters chosen is bandwidth efficient, so that it does not require narrow 
filtering and consequently offers constant or near-constant envelope. 


4 Organization of this report 

This Final Report is organized as follows. Chapter 1 provides the motiva- 
tion for this work, as well as an overview of the channel model for which 
our system has been analyzed zmd designed. Chapter 2 reviews continuous- 
phase modulation, with emphasis on the spectral properties of continuously 
phase-modulated signals. Chapter 3 reviews trellis-coded modulation, and 
discusses the interactions between TCM and CPM, as well as the possible 
advantages resulting from it. Chapter 4 deals with suboptimum detection 
of trellis-encoded CPM. In particular, two extreme cases of suboptimum 
detection, viz., coherent and noncoherent, are analyzed. Differential de- 
tection of trellis-encoded CPM is the subject of Chapter 5. Chapters 6 
and 7 describe the design and the analysis of two circuits that are needed 
for the proper operation of a differentially-detected, trellis-encoded CPM 
transmission system. In particular. Chapter 6 considers the effects of a 
Doppler frequency shift, and its removal. Chapter 7 considers a circuit for 
the recovery of the timing information from the received signal. SimiJation 
results are included in Chapter 8, while Chapter 9 contains the conclusions 
drawn from the preceding body of work. Two appendices describe the sim- 
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ulation package written during the work and the hardware implementation 
of the transmitter and receiver^. 


^We observe here that this Final Report is not the result of a simple superposition of 
our Quarterly Reports. For legibility's sake, some discussions that are less relevant to our 
presentations, as well as a number of charts showing results of lesser importance, have 
been omitted, and the whole material has been reshuffled. 
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Chapter 1 

Introduction 






In this Report we consider a satellite-based mobile communication scheme 
bcised on continuous- phase modulated (CPM) signals used in conjunction 
with trellis- coded modulation (TCM). 

In satellite-based land mobile communication systems both bandwidth 
and power are limited resources. In fact, these systems employ frequency- 
division multiple access (FDMA) with a given channel spacing (say, 5 KHz), 
and the fraction of out-of-band power should be very small to prevent inter- 
ferences to adjacent channels. On the other hand, the satelhte distance from 
earth, its power limitations, and the need for low-cost (and hence low-gain) 
mobile antennas put a serious limit on the power resources. Additionally, 
the fading environment of mobile communication further limits the power 
efficiency of the system. 

In a bandwidth- and power-limited environment, a bandwidth- and power- 
efficient coding/modulation scheme must be used. Trellis coded modulation 
(TCM) offers an attractive scheme. It combines the choice of a higher-order 
modulation scheme with that of a convolutional code, while the receiver, 
instead of performing demodulation and decoding in two separate steps, 
combines the two operations into one. As a result, the reliabihty of a digital 
transmission system is increased without increasing the transmitted power 
nor the required bandwidth. By using TCM, simple schemes can be de- 
signed that achieve significemt power gains (from 3 to 6 dB) without any 
bandwidth expansion (see, e.g., [60,80].) 

Due to the strictly bandlimited environment created by the mobile satel- 
lite channel, the signals to be used in conjunction with trellis codes must be 
chosen carefully. Besides having excellent spectral characteristics, the sig- 
mls used should have constant envelope if nonlinear amplifiers are used for 
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better power efficiency. A class of bandwidth-efficient signals that satisfies 
both constraints is offered by continuous-phase modulated (CPM) signals, 
based on phase modulation where phase continuity is introduced to reduce 
the beindwidth occupancy. 

The synergy between TCM, which improves error probability, and CPM 
signals, which provide constant envelope and low spectral occupancy, is ex- 
pected to provide a satisfactory solution to the problem of transmitting on 
mobile satellite channels. 

In principle, two implementations of this idea are feasible. The first one 
takes advantage of both the bandwidth efficiency and the power efficiency 
of CPM codes, by using a receiver which combines the trellis structure of 
TCM and that of CPM. In this situation, TCM and CPM can be integrated 
in a single entity (see [51] and the references therein), but the number of 
states necessary for a trellis representation of these signals is given by the 
number of TCM states times the number of CPM states. As this number 
can grow very large, the complexity of the receiver may become quickly im- 
manageable, and a suboptimum solution should be devised. We can trade a 
decrease in complexity for a decrease in power efficiency (but not in band- 
width efficiency) by giving up maximum-likelihood decoding of the CPM 
signals, which are instead demodulated symbol-by-symbol. By doing this, 
the power efficiency of CPM codes is not exploited: we only take advsintage 
of their spectral properties. 

Hereafter we consider three solutions: coherent, differentially-coherent, 
and non-coherent symbol-by-symbol detection. Since differentially coherent 
detection appeared to us to be the most promising technique, most of this 
Report is devoted to a complete description of a CPM/TCM system based 
on this kind of detection. 


1.1 The channel model 

Besides additive Gaussian noise, which is the standard environment for the 
analysis of coding schemes for the transmission of digital data or speech, 
there is a number of additional sources of performance degradation that 
must be taken into accoimt to assess the merits of a proposed transmission 
scheme for mobile satellite channels. The most important among them are 
[61,64]: 

• Doppler shifts. They are due to mobile vehicle motion. If differenticd 
detection is used, the information-bearing phase turns out to be shifted 
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Velocity 

(MPH) 

Frequency 
shift (Hz) 

Phase shift 
(degrees) 

30.0 

67.04 

10.06 

60.0 

134.08 

20.11 

90.0 

201.13 

30.17 

120.0 

268.17 

40.22 


Table 1.1: Doppler frequency and phase shifts. 

by an amount 27 r/dT,, where l/T, is the data symbol rate and fd is 
the Doppler frequency shift, which for operation at L-band can be 
expected to be up to 200 Hz* At a rate of 2400 symbols per second 
the corresponding phase shift is Table 1.1 shows some values 

of the Doppler shift obtained for various velocities, carrier frequency 
1.5 GHz, and symbol interval T, = 1/2400. 

• Fading and shadowing* The transmitted radio signal reaches the 
receiver through different paths caused by reflections from obstacles, 
yielding a signal whose components, having different phases 2 md ampli- 
tudes, may either reinforce or cancel each other. Shadowing is caused 
by the obstruction of radio waves by buildings, trees, and hills. 

• Adjacent channel interference* The 5-KHz mobile channel used 
for transmission operates in a channelized environment. As a result, 
signals suffer from interference from signals occupying adjacent chan- 
nels. 

• Channel nonlinearities* Primarily because of the high-power am- 
plifier in the transmitter, operated at or near saturation for better 
power efficiency, the channel is inherently nonlinecir. 

• Finite interleaving depth* In order to break up the error bursts 
caused by amplitudes fades of duration greater than symbol time, en- 
coded symbols should be interleaved. Now, infinite interleaving pro- 
vides a memoryless channel, but in practice the interleaving frame 
must be limited. In fact, for speech transmission the total coding/de- 
coding delay must be kept below 60 ms in order not to cause per- 
ceptual amnoyance. If the depth of interleaving cannot be larger than 
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the m^Yirmim fade duration anticipated, this causes 
degradation. 


a performance 


Chapter 2 


Continuous Phase 
Modulation (CPM) 


In this Chapter we provide an overview of the basic features of continuous- 
phase modulated signals. Additional details can be found in [55,58]. 

A continuous-phcise modulated signal is defined by 


s(t,a) 



{2rfot + 


( 2 . 1 ) 


where E, is the symbol energy, T, is the symbol time, and /o the carrier 
frequency. The transmitted information is contained in the phase 


OO 

6{t,a) = 2irh ^ a„«(t - nT,) (2.2) 

n=— OO 


with q{t) the phcise- shaping pulse given by 

?(0 = / (2-3) 

J — OO 

g{t) is the frequency pulse with finite duration, i.e., ^(t) is nonzero only for 
0 < t < LTsy L the pulse length. The value of L contributes to the taxonomy 
of CPM signals, namely: 

• L = 1: Pull- response signals 

• L > 1: Partial- response signals. 
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It is common to assume 


r g{r)dr= 1/2. 

J — 00 


In (2.2), 

denotes the symbol sequence sent to the CPM modulator. The symbols an 
take values ±1, ±3, • • ■ , ±{M - 1), where M = 2"*, m a positive integer. The 
parameter h is called the modulation index, and we shall assume 

i 

with p,q relatively prime integers. Notice that the change in the instanta- 
neous frequency, i.e., 

fo + 9(t, a)/2x 

is proportional to the modulation index. The maximum phase change over 
a symbol interval is Trh[M — !)• 

When the duration of the frequency pulse ^(t) is greater than 1 (partial 
response CPM), the phase function d(t,a) during the symbol interval may 
be written in the form 

e{t, a) = 2xh ^ ajq{t - jT,) + 9n, nTs<t<{n + 1)T, (2.4) 

j=n-L+l 

where 

n^L 

^ aj mod 2x. (2*5) 

i=-oo 

It can be seen from (2.5) that can take q different values, namely 

Thus, we can say that the phase function during any given interval depends 
on the actual tr£insmitted symbol and on the state of the modulator, 
defined as the vsilue taken by the vector 

(an-l»fln-2»”*»“n-L+l>^n) (2*7) 

For each state there axe M signal trajectories, and the total number of 
distinct signal paths over any T^-second interval is qM^, The total number 


CHAPTER 2. CONTINUOUS PHASE MODULATION (CPM) 


16 


of states is As a special case, for ftiU- response signaling there are 

Mq signal paths, and q states: this reduction in the number of paths and 
states, and hence in the complexity of the modulator-demodiilator pair, is 
traded off for an inferior spectrum. 

2.1 Optimum detection of CPM Signals 

Optimum (maximum-lilcelihood sequence) estimation of CPM signals in- 
volves maximization of the probability density function for the observed 
signal conditioned on the symbol sequence a. Under the assumption that 
the only distuirbance adfecting the received signal is an additive white Gaus- 
sian noise process n(t), i.e., that r(f) = s(t) -t- u(t), optimum detection is 
equivalent to meiximizing the quantity 



J(a) = 

too 

1 r{t)s{t,a)dt. 

J — OO 


If we define 

•fn(a) = J 

(n-l-l)T. 

r(t)s(f,a)dt, 

•OO 

(2.8) 

we can write 

^(a) = 

<^n-l(a) + ^n(a). 

(2.9) 

where 

/■(n-M)T, 

2n(a) = 

J nTs 

r(t) cos[2ir/ot -|- 0(t,a)]dt, 

(2.10) 


CcJled the branch metric, is the correlation between the received signal and a 
reconstructed version of the transmitted signal over the n-th symbol interval. 
By using (2.8)-(2.10) it is possible to compute J(a) recursively. The Viterbi 
algorithm chooses the sequence a that maximizes Jn(a). 

Because of (2.4) we can zJso write 

y*(n+l)Tj 

Zn{a,0n)= / r(f)cos[2x/ot -b d(t,a„) + (2-H) 

JnTs 

where 

n 

d{t,&n) = 2irh ^ ajq{t - jT,). 

j=n-L+l 

The receiver computes Z„(an,^n) for all possible sequences 


an = (®n) ■ • • jOn-L-n) 
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and aU g possible This shows that can take on gM possible v^ues. 
The computation of these values requires feeding the observed signal r(t) 
into a bank of qM^ matched filters which correlate the received signal over 

one symbol interval. 

A maximum-likelihood detector provides a bit error probabihty that for 
high signal-to-noise ratios can be roughly approximated by 



where d? is the minimum integral-squared difference between any two sig- 
nals corresponding to phase trajectories that are merged till a certam time, 
then split and remerged at an arbitrary time later. 

2.2 Computation of Power Spectrum 

We consider here the evaluation of the power spectral density of CPM sig- 
nals. The method used is based on a technique developed by Auhn an 
Sundberg [56], and computing the spectrum by Fourier transformation o 

the autocorrelation function of the modulated signal. . ,. 

The essence of the method is the following. First, the CPM si^al is 
multiplied by its shifted version and its expectation is taken to obtam the 
autocorrelation. As the CPM signal process is not wide-sense stationary 
(but rather cyclostationary), a time-averaged version of the autoco^elation 
function is determined, then its Fourier transform is taken to get the power 
density spectrum. 

First, let us rewrite (2.1) in a complex form as follows: 

,(t,a) = (2.13) 

where F = E,ITs is the average signal power, and 71 denotes “real part . 
The other symbols are the same as in (2.1). We assume here, for simpUcity, 
that the data symbols are independent: this assumption is not consistent 
with our model if trellis-encoded CPM is considered, but we assume that 
the effect of TCM on the power spectrum is negUgible. 

Let M denote time average (over the interval (0,T,)), ^d E denote 
expectation with respect to the random variables a„ representing the trans- 
mitted data. The (stationary) autocorrelation function of the CPM signal 
is then, under the assumption foT, > 1, 

r(r) = ME{s(t -I- r,a)a*(f,a)} 

= P7^{i^(T)e^2x/orJ^ 
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Here, 


i?(r) = 

_ M £^gi2irAai[q(t+T-«T,)-9(t-iT,)]j._ 


Now, 

E{gi2?r/io,[9(t+T-iT,)-g(t-tT,)]} 

M/2-1 

_ ^ J2nh(2k+l)[q{t+T-iT,)-q(t-iT,)] 

~M^ 

k=-M/2 

Moreover, noting that q{t) = 1/2 for t > XT,, and q(t) = 0 for t < 0, and 
letting 

T = "1“ TTiTg , 0 ^ T ^ rn — 0, 1, 2, • • • , 

it can be shown that 


R{t) = R{T' + mT,) 

rT, vi+i 


Mf 2-1 


1 rT, 1 ‘ 

= yL n s e 

JO i=i_L fc=-A//2 


SO that the power spectrmn takes the expression 


h(2k+l)[q{t-^r^-{i^m)Ts)-q(t-iTs)]^^ 


Oif) 

where 




R{r)e- 




rTs 


M 


1 - Jo 

R{t + 

(2.14) 

Af/2-1 

e-'’"**- 

k=-M/2 

(2.15) 


To compute R{t) effectively and accurately, we subdivide the integration 
interval (0,Tj) into n subintervals, and use a 5-point Gaussian quadrature 
formula to compute the integral over each subintervad. The number n will 
be chosen to be large enough so as the limitation in accuracy comes from 
computational roundoff errors rather than from numerical approximations. 
After computing the autocorrelation function, the fast-Fourier transform 
algorithm is used to obtciin the power density spectrum. 

Power density spectra, in the form of fractional out-of-band power (in 
dB) were calculated for several different pulse shapes. The results for four 


CHAPTER 2. CONTINUOUS PHASE MODULATION (CPM) 19 

of them (namely, LREC, LRC, HCS, and CRC for L = 1,2,3, and 4) are 
shown in the next figures. SpecificaUy, Fig. 2.1 shows the expressions for 
the frequency pulses in the four cases considered here. Fig. 2.2 shows the 
corresponding expressions for phase pulses. Frequency pulses and phase 
pulses are tabulated in Fig. 2.3 and Fig 2.4, respectively. Figs 2.5 to 2.12 
show the fractional out-of-beind powers. 
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Figure 2.4: Phase pulse g{t) for LREC, LRC, HCS, and CRC. 
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Figure 2.5: Fractional out-of-band power (in dB) of octonary CPM signals 
with modulation index h = 1/8 for LREC pulse shape, L — 1, 2,3,4. 
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Normalized Frequency 


Figure 2.6: Fractional out-of-band power (in dB) of octonary CPM signals 
with modulation index h = 1/8 for LRC pulse shape, L — 1,2,3, 4. 
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fT, Normalized Frequency 


Figure 2.7: Fractional out-of-band power (in dB) of octonary CPM signals 
with modulation index h = 1/8 for HCS pidse shape, L = 1,2, 3,4. 
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Figure 2.9: Compcirison of fractional out-of-band power (in dB) of octonary 
CPM signals with modulation index h = 1/8 for different pulse shapes, 
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Fieuie 2.11: Comparison of fractional out-of-band pow^(in dB) of octonary 

* 1 1 h 1 /R for different pulse stapes, 

CPM signals with modulation index — / 
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Figttre 2.12: Compajison of fractional out-of-bfind power (in dB) of octonary 
CPM signals with modulation index h = 1/8 for different pulse shapes, 
L = i. 




Chapter 3 

Trellis Coded Modulation 
(TCM) 


In this Chapter, we present a review of the main features of trellis-coded 
modulation (TCM) as used for digital data communication with the purpose 
of gaining noise immunity over uncoded transmission without altering the 
data rate. 

For the purpose of this discussion, we shall cissume a discrete-time, 
continuous-amplitude model for the transmission of data on the additive 
white Gaussian noise channel. In this communication model, the messages 
to be delivered to the user are represented by vectors in an IV- dimensional 
Euclidean space R.^, called the signal space. When the vector x is trans- 
mitted, the received signal is represented by the vector 


z = X -1- y. 


where y is a noise vector whose components are independent Gaussian ran- 
dom variables with me 2 m zero aind the same variance Nq/2. The vector x is 
chosen from a set Q consisting of M signal vectors and which will be referred 
to as the signal constellation. The average square length 




xefi 


will be referred to eis the average signal energy. 

Consider now the transmission of a sequence of K signal vec- 

tors, where the subscript i denotes discrete time. The receiver which min- 
imizes the average error probability over the sequence operates as follows. 


32 


CHAPTER 3. TRELLIS CODED MODULATION (TCM) 


33 


It first observes the received sequence Zq, • • • , then it decides that 

Xo, . . . was treinsmitted if the squared Euclidean distance 

i=0 

is minimized for x; = X,, i = 0, — 1, In words, the sequence 
Xo, is closest to the received sequence thcin any other zdlowable 
signal vector sequence. The resulting sequence error probability, as well as 
the symbol error probability, are upper bounded, at least for high signal-to- 
noise ratios, by a decreasing function of the ratio <^a/No, where is the 
minimum squared Euclidean distance between two allowable signal vector 
sequences. 

Uncoded Transmission 

An important special case occurs when the signal vectors form an inde- 
pendent sequence. In this case, the allowable vector sequences are all the 
elements of and hence is minimized by minimizing separately the 
individual terms of (3.1), i.e., |Z,- — Xjp for x^ G and i = 0,. . . yK — 1. 
The performance of this “symbol-by-symbol receiver” will then depend on 
the minimum distamce 


dLn = |x' - x' 


«|2 


x'#x‘ 


as x^,x" run through fl. In fact, the symbol error probability is upper 
bounded (and for high signal- to-noise ratios well approximated) by 


P(e)<M^»fc(^) 


(3.2) 


(see, e.g., [58, page 152]). 

With this model, the problem of designing an efficient communication 
system is that of choosing a set of vector signads such that the minimum 
distance between any two of them is a maximum, once the quantities M,N, 
and £ ^l^e given. It is convenient to define two quantities that are useful in 
comparing different constellations, namely their information rate (measured 
in bits /dimension): 

logM 


(3.3) 
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and their normalized squared minimum distance 

^2 = ^logM. (3.4) 

The first parameter is also referred to as the “bandwidth efficiency”, because 
it represents the ratio between the number of information bits carried by a 
single signal in the constellation and the number of dimensions, which is 
roughly proportional to the transmission bandwidth- The second parameter 
can be referred to as the “energy efficiency” of the signal constellation: in 
fact, observe that the upper bound (3.2) can be rewritten in the form 

where ^ 

^ logM 

represents the average energy per bit. It is seen from (3.5) that the same 
P{e) can be achieved with a a smaller signal-to-noise ratio if # is 

larger. 


3.1 The Concept of TCM 


TCM is a coded system: this means that it tries to achieve an improvement of 
system performance by introducing a redundancy, i.e., an interdependency 
between signal vectors. This is equivalent to restricting the transmitted 
sequences to a subset of Now, if we do this, the transmission rate 

will also be reduced. To avoid this unwanted reduction, we may choose 
to increase the size of fl. For ex£unple, if we change 0 into D and 
M into M' > M, and we select sequences as a subset of fl' , we 
can have sequences which are less tightly packed, and hence increase the 
minimum distance between them. In conclusion, we get a minimum distance 
dfree between two possible sequences which turns out to be greater than the 
minimum distance between signals in fl, the constellation from which 
they were drawn. Use of maximum-Ukelihood sequence detection wiU yield 
a “coding gain” <^i^/d^n 5 the additional power necessary for signaling 
using the alphabet 0' instead of fl. We define the energy gain of a TCM 


scheme as 
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where £ and £' denote the average energies spent to transmit with uncoded 
and coded tramsmission, respectively. 

The introduction of interdependencies between signals as a way of in- 
troducing the wanted redundancy into the sequence set is one of the basic 
ideas underlying trellis-coded modulation. (Another one is set partitioning, 
which win be described later). 

We assume that the signal x„, transmitted at discrete time n, depends 
not only on the source symbol Un transmitted at the same time instant (as 
it would be with memoryless modulation), but also on a fmite number of 
previous source symbols: 

= /(®n,an-l)-’->®n-L)- (3-6) 


By defining 


— (®n— 1 1 " ’ * 1 l) 


(3.7) 


as the state of the encoder at time n, we can rewrite (3.6) in the more 
compact form 


3Cn — /(<*n» <^n) 
O’n+l = <7(<*n,<^n)* 


(3.8) 

(3.9) 


The last two equations can be interpreted as follows. The fimction /(•, •) 
describes the fact that each channel signal depends not only on the corre- 
sponding source symbol but also on the parameter <Tn. In other words, at 
any time instant the transmitted signal is chosen from a constellation which 
is selected by the value of <r„. The function ?(•, •) describes the memory part 
of the encoder and shows the evolution of the modulator states (see Fig. 3.1). 
Here we shall Msume that the functions / and g are tune-invariant, although 
it is possible to consider time- varying TCM schemes as well. 


3.1.1 Trellis Representation 

For a graphical representation of the frmctions / and g it is convenient to use 
a trellis. The values that can be taken by <r„, the encoder state at time n, 
are the nodes of the trellis. With each source symbol we associate a branch 
which stems from each modulator state at time n, and reaches the encoder 
state at time n -|- 1. The branch is labeled by the corresponding value of /. 
The trellis structure is determined by the function g, while / describes how 
channel signals are associated with each branch along the trellis. 
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Figure 3.1: A general model for TCM 
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If the source symbols axe ilf-ary, each node must have M branches stem- 
ming from it (one per source symbol). This implies that in some cases two 
or more breinches connect the szime pair of nodes: when this occurs, we say 
that parallel transitions take place. 

Fig.3.2 shows am example of this representation. It is assumed that the 
encoder has four states, the source emits binary symbols, and a constellation 
with four signals 0,1, 2, 3 is used. The distance properties of a TCM scheme 
can be studied through its treUis diagram. 


3.2 Some Examples of TCM Schemes 

Here we describe a few examples of TCM schemes based on their trellis 
diagrams. We first consider the transmission of quatemairy source symbols, 
i.e., 2 bits per symbol. With uncoded transmission a channel constellation 
with M = A would be adequate. We shall examine TCM schemes with 
M' = 8. 

Let us assume M-ary coherent PSK transmission. With M = 4 we get 



2 , 


a figure which will be used as a reference to compute the coding gain of PSK- 
based TCM. We use TCM schemes based on two quaternary constellations 
{0, 2,4,6} and {1,3, 5,7} shown in Fig. 3.3. We have 


€' = 


d'2 

4sin2f 


Consider first a scheme with two states as shown in Fig. 3.4. If the encoder is 
in state 5i, constellation {0, 2,4, 6} is used. If it is in state 52, constellation 
{1, 3,5, 7} is used instead. (Notice the presence of parallel transitions). 

Let us compute the free distance of this scheme. This can be done by 
choosing the smallest between the distance among signals associated with 
parallel transitions, and the distance associated with a pair of paths in the 
trellis that originate from a common node jmd merge into a single node at 
a later time. The pair of paths giving the free distance is shown in Fig. 3.4, 
and we have, by denoting d(i, j) the Euclidean distance between signals i 
and j: 

^ = l[d2(0, 2) + d^(0, 1)1 = 2 + 4 sin^ ^ = 2.586. 
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Figure 3.2: An example of trellis describing a TCM scheme with four states 
and four signals used to transmit from a binary source. 
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Hence, we get a coding gain 


7 = 


2.586 

2 


= 1.293 ={► 1.1 dB. 


Let us now use a TCM scheme with a more complex structure to increase 
the coding gain. With the same constellation of Fig. 3.3, take a trellis with 
4 states as in Fig. 3.5. We associate the constellation {0,2, 4, 6} to states 
Si and S 3 , and {1, 3,5, 7} to 52 and S^. The pair of paths giving has 
length 1 (a parallel transition) and is shown in Fig. 3.5. We get 

^ = d2(0,4) = 4. 


and hence 

7 = — = 2 => 3 dB. 

^ 2 


A fiirther step can be taken by choosing a trellis with 8 states as shown in 
Fig. 3.6. For simplicity, the four symbols associated with each node are used 
as node labels. The pair of paths leading to dfre* is also shown and yields 



■i [(P (0, 6) + d^(0, 7) + d^(0, 6)] = 2 + 4 sin^ I + 2 = 4.586, 

O 


and hence 


7 = 


4.856 

2 


= 2.293 


3.6 dB. 


Consider now the transmission of 3 bits per symbol and AM-PM schemes. 
The octonziry constellation 


{0,2,5, 7,9,10,13, 15} 


of Fig. 3.7 will be used £is a reference uncoded scheme. It yields 

(P. 

= 0 . 8 . 

A TCM scheme with 8 states and based on this alphabet is shown in Fig. 3.8. 
The sub constellations used sue 


{0,2,5, 7,9,10,13, 15} 


and 


{1,3,4,6,9,11,12, 14}. 
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Figure 3 


.6: A TCM scheme based on an 8-state trellis, M = 4, and M' — 8 
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0,2,5,7,8,10,13,15 

I, 3,4,6,9,11,12,14 

5.7.0. 2.13.15.8.10 

4,6,1,3,12,14,9,11 

10.8.15.13.2.0. 7.5 

II, 9,14,12,3,1,6,4 

15.13.10.8.7.5.2.0 

12.14.9.11.6.4.3.1 



Figure 3.8: A TCM scheme based on an 8-state trellis and AM-PM constel- 
lations, with M = 8 cind M' = 16 
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£’ = 2.5 d'^. 

I[d2(10,13) + d'*(0,l) + <i"(0,5)] 

1[0.85' + 0.45' + 0.85'] 
o 
2 

and hence 

7 = — = 2.5 => 3.98 dB. 

' 0.8 

Set Partitioning 

Consider the determination of dfree- This is the distance between the signals 
associated with a pair of paths that originate from an initial split, and, after 
L time instants, merge into a single node as shown in Fig. 3.9. Assume first 
that the free distance is determined by parallel transitions, i.e., i = 1. If A 
denotes the set of signals associated with a given node, then equals the 
minimum distance among signals in A. 

Consider then I > 1. With A,B,C,D denoting subset of signals as- 
sociated with each branch, and d(JT, Y^) denoting the minimum Fuclidean 
distance between one signal in X and one in Y , will have the expression 

+ * •• + d^(C, I?). 

This implies that, in a good code, the subsets assigned to the same origi- 
nating state or to the same terminating state (“adjacent transitions”) must 
have the largest possible distance. To implement these rules, the technique 
suggested by Ungerboeck [80] and called set partitioning is the following. 
Set partitioning has been described as “the key that cracked the problem of 
constructing efficient coded modulation techniques for band-limited chan- 
nels.” 

The constellation is successively peu*titioned into 2, 4, 8, , sub- 

sets with size M'/2, M74, M'/S, . . . , having progressively larger minimum 
distances • (see Fig.3.10). Then, 

1. Parallel transitions are assigned members of the same partition. 

2. Adjacent transitions are assigned members of the next larger partition. 
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A 



B 

L = 1 



L> 1 


Figure 3.9: A pedr of splitting and remerging paths for 
transitions) and L > 1. 


<^n + L 


= 1 (parallel 
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3.3 Combining TCM With CPM 

If CPM sigii£ds are combined with an external convolutional encoder, or 
equivalently they are used as the signal constellation to be used in a TCM 
scheme, a further improvement can be obtained. This new scheme is ob- 
tained by observing that at the output of the trellis encoder we get a mtilti- 
level signal which in turn can be used as the input to the continuous-phase 
modulator. The design of the coding scheme and of the modulator scheme 
should be performed jointly in order to maximize the Euclidean distance 
resulting from the combination of the two. As observed in [76], to maximize 
the Euclidean distance of the coding/modulation scheme, the trellis encoder 
should reduce the connectivity of its trellis in such a way that: 

• The Euclidean disteince between signals leaving the same state is max- 
imized. 

• The Euclidean disteince merging into the same state is maximized. 

To implement TCM/CPM, we may want to take advantage of both the 
bandwidth efficiency and the power efficiency of CPM codes. In this case, the 
receiver combines the trellis structure of TCM with that of CPM, and TCM 
and CPM can be integrated in a single entity (see [51] and the references 
therein). The number of states necessary for a trellis representation of these 
signads, and hence for their demodulation, is the product of the number of 
states needed by TCM and the number of states needed by CPM. Since this 
number can grow very large, a suboptimum solution should be sought, which 
trades a decrease in complexity for a decrease in power efficiency (but not in 
bandwidth efficiency). This is obtained by using a differential demodxilator 
or a discriminator instead of a coherent receiver. 


Chapter 4 

Coherent vs. noncoherent 
CPM detection 


To keep the demodulation problem separate from the decoding problem, 
we consider estimating the CPM phases symbol-by-symbol, and using these 
estimates to build up a metric for the decoder. Provided that the number 
of decoder states is finite and the metric chosen is additive, we can use a 
Viterbi algorithm for demodulation, the number of states being that related 
to the TCM scheme only. This procedure obviously entails a loss of opti- 
mality, which is traded against a manageable receiver structure. However, 
as mentioned before, the spectral properties of CPM are preserved. 

In this Chapter we consider coherent and noncoherent symbol-by-s 3 ncnbol 
demodulation of CPM signals. Coherence can be achieved by taking advan- 
tage of the synchronization properties of CPM [38,62], sJthough, as observed 
in [57,52], as the spectrum of the transmitted sign 2 il becomes narrower, the 
problem of finding the exact phase of the carrier increases. To avoid the 
problem of finding the phase of the carrier, one viable solution is the use of 
noncoherent detection. 

Suboptimum coherent and noncoherent detection of CPM (in the special 
case of CPFSK) was considered in [48] for uncoded transmission. In this 
Chapter, consideration of suboptimum schemes involves the derivation of a 
metric to be used with the Viterbi decoder for TCM. Thus, the appropriate 
performance measure in our case is the cutoff rate Rq [47] of the channel 
generated by CPM in conjunction with its detection scheme. The use of 
Rq to demonstrate the efficiency of a coding scheme in the case of a fading 
channel, typicad of mobile satellite and terrestrial mobile radio applications. 
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was advocated in [45]. 

The orgainization of this chapter is the following. We first summarize our 
signal and channel model. Symbol-by-symbol coherent detection of CPM 
is first examined, then incoherent detection is covered in Section 4. The 
Appendix reviews the basic steps involved in the computation of Rq for 
fading channels. 


4.1 Signal and channel models 

The communication scheme considered in this section is shown in Fig. 4.1. 
The source data 

u = uo,ui,--- 

are first sent into an encoder component, whose output is the symbol se- 
quence 

a = ao,ai,**< 

The encoder includes a conventional linear convolutional encoder with rate 
m/(m -I- 1), followed by an M-ary mapper. The convolutionally encoded 
output sequence consists of binary data symbols. The output of the mapper 
is the sequence a of M-ary symbols, generated according to some specific 
mapping rule. We assume that M is a power of 2, and that the symbols an 
take values ±1,±3, •- • j±(M — 1). 

The coded sequence a is first interleaved (by using a block or convolu- 
tional interleaver [41, pp. 347-349]), to produce a'. This new sequence is 
then sent to the continuous-phase modulator, which outputs a signal whose 
complex envelope is 

s{t, a', <f>o) = ^^eip j[e{t, a') + ^o] (4.1) 

where 4>o is the carrier phase, E, is the signal energy per symbol, and T, is 
the symbol time. The energy per bit is Ei, — F./logj M. The treinsmitted 
information is contained in the phase 


9{t, a') = 2irh a',q{t - kT,) (4.2) 

fc =0 

with q{t) the phase-shaping pulse given by q{t) = /J g{T)dT, and g{t) is the 
frequency pulse with finite duration, i.e., g{t) is nonzero only for 0 < t < 
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LTs, L the pulse length. The parameter h is called the modulation index, 
and we assume it to take on a rational value. The maximum phase change 
over a symbol interval is vh{M — 1). 

A CPM scheme is defined by selecting the values of M, h, and the fre- 
quency pulse shape. Some of the most popular pulse shapes are LREC (for 
L = \ often referred to as Continuous- Phase Frequency Shift Keying, or 
CPFSK), Tamed Frequency Modulation (TFM), Raised Cosine (LRC), Half 
Cycle Sine (HSC), Convolved Raised Cosine (LCRC), Gaussian Minimum 
Shift Keying (GMSK), Spectrally Raised Cosine (LSRC). (For further de- 
tails, see [55,58].) In this Chapter we consider Full Response (i.e., L = \) 
CPM. 

In our analysis, we assume that the channel includes multipath reflec- 
tions in addition to Gaussian noise. Thus, if r(f) and i/(t) denote the signal 
observed at the channel output and a complex Gaussian noise process, re- 
spectively, we have 

r(t) = p(t)s(<,a',</>o) + »'(0 

where it is assumed that r(t) is observed during the time interval (0, acTj), 
and fi(t) represents the (normalized) random fading amplitude. The prob- 
ability density function of p(t) depends on the model chosen for the fading 
process: we have 

p(p) = 2pe~^\ p>0 

for the Rayleigh model, 

p{p) = 2p{l + -1- K)), p>0 

for the Rician model In the latter expression, K denotes the ratio of 
power in the line-of-sight component to that in the difftise or noncoher- 
ent component, md 7o(-) the modified Bessel function of order 0. Note that 
E[ipp] = 1 for the Rayleigh and (unshadowed) Ricim model. Note also that 
the Rayleigh fading model is the limiting case of the Rician model when 
K = 0. 

4.1.1 Computation of Rq 

Here we summarize the relevant steps in the derivation of the cutoff rate 
Rq of the memoryless chjinnel generated by CPM in conjunction with an 
interleaver /deinterleaver pair. Consider the transmitted codeword x;v with 
block length N, Let us assume that decoding is based on the observation of 
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the received iV- tuple r;v &nd on the use of the additive metric 


N 


m(rN,xiv) = ^rn(ri,Xi), 


( 4 . 3 ) 


1=1 


where r;, xi are the components of riv, xu, respectively. The decoding rule 
is then: choose X;v if 

m(r\r,iiv) = max m(r/\f,XN}- 

In this situation, the pairwise error probability, i.e., the probability that 
xjv have a larger metric than the transmitted x^v, can be bounded 

from above by using Chemoff bound. Under our assumption that the chan- 
nel is memoryless, we have 


P[x^r -*■ Xiv] < E 
N 


I E.'Ii -^•)] |xiv I 


Now, define the quantity Dx as follows: 

g-r>x(2:.x) ^ E I , 


( 4 . 4 ) 


so that ^ 

p[xn^±n] < 

1=1 

Next, apply a random coding argument. With the assumption that 
the code alphabet consists of M distinct signals, randomly select the 2N 
signals forming xyv and xjv^ where each signal is independently selected 
with uniform probability. The average pairwise error bound is then 


P[X/^ X;V] < 


N 

JJ ^-Dx(xiyX^) 

i=l 

N 

JJ Q-Dx{xi,xt) 
1=1 

2-Ro{\)N 
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where we have defined 

Ro{\) = - log 2 (4.5) 

and the random coding average has the explicit expression 

g-Dx(x,x) — g~D\{x,x)^ ^4 5 ) 

X X 

The value of A should be chosen in order to minimize the Chemoff bound, 
i.e., to maximize Rq{X). If this is not feasible, a suboptimum bound will be 
obtciined. 


Fading channel 

With a channel affected only by additive complex white Gaussian noise n 
with two-sided power spectrzil density 21Vo) fading (represented by a 
random variable p), if perfect state information is assumed we may use the 
metric 

m(r,a:) = - || r - p* |p, 

where r = px + n. In this case, by observing that for a complex Gaussian 
random variable n with mesm zero, variance 2iVb and independent real and 
imaginary parts we have 


a simple calculation shows that, by choosing for A its optimum value l/4JVo, 




= E |< 





where the expectation is taken with respect to the random variable p repre- 
senting the fading effect. We get the expression 






)■ 


(4.7) 


4.2 Suboptimum coherent detection 

The detection problem is as follows. The optimum (maximum likelihood 
sequence) detector must find the sequence of source symbols u which is 
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mapped into the sequence a' that maximizes the likelihood ftmction 
A[r(t)la',p(t),<;6o] = |r(t) - 

= exp I II r(t) - p(t)s(t, a\^o) ||^ dtj . 

Now, let us divide the total observation interv 2 d into k subintervals of dura- 
tion Ty seconds each. Then we can write 

1 

II r(t) - p(t)s(t, a\<f>o) |p= XI II “ Pk^(°ki<f>k7^) 11^ • 

fc=o 


The time dependence has been omitted for simplicity, and pk denote r(t) 
and p(t) in the interval (ibT„ (ib + 1)T,), and, in the same time interval, 


where 

it-i 

<pk “ ^ ^ k — 1, . . . , /C, 

n=0 

is the value of the phase accumulated up to the fc-th time interval. The 
latter is also cedled the “phase state” of the CPM signal. 

If we assume that the receiver achieves perfect coherence, the value of 
is known, and hence, without loss of generality, we shall assign it the value 0 
and exclude it from consideration as well as from our notations. Disregard, 
for a moment, the interleaver/ deinterleaver pair. In this situation, optimum 
coherent detection of the received sequence would be obtained by using in 
the TCM detector the additive metric 

m[r(«),3(i,a')] = X II '’fe “ Pfc^(4><Afc) |P • 

fc=o 

Tti the Viterbi-adgorithm parlance this corresponds to associating the branch 
metric || 11^ trellis brainches corresponding to the 

transmitted symbol a'^.. Now, it is apparent that the metric values depend 
on 4>k as well, i.e., on the state of the CPM modulator when the signal 
s{a'k,<i>k) was transmitted. To take this fact into account, we should expand 
the TCM trellis. Let us denote by 5 tcm the number of TCM states and 
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by 5cpm the number of CPM states. The latter is equal to the number of 
distinct vsdues (mod 2ir) of the accumulated phases, i.e., if h = pjq is the 
modulation index, with p emd q relatively prime integers, equal to either 
9 (if p is even) or 2q (if p is odd) [58, p. 169]. Then every state in the 
original TCM trellis must be transformed into 5cpm states, one for each 
value of 4>k- The restilting “product trellis” has 5 tcm^CI^ states, and this 
is the number that determines the computational complexity of the “single- 
approach” Viterbi algorithm. 

To avoid this state growth and to be able to perform interleaving/ de- 
interleaving (which is feasible only if we separate the operations of CPM 
demodulator and TCM decoder), we should get rid of the term <f>k, possibly 
by using non-coherent or differential detection (see, e.g., [36]). Here we 
consider a coherent scheme, based on the following idea. The CPM signal is 
sent to an optimum coherent detector, based on the Viterbi algorithm and 
which outputs only the sequence of states, i.e., the values of the sequence 
<f>k- These phases aire then removed from the received signal, which turns 
out to be, in the interval (fcT„ {k + 1)T,), 

/or 

r{t) = J -jr - kT^) + u{t). 

After deinterleaving, the original order of the sequence a is restored, and 
hence the metrics 1| Vk - Pk^{(^k) |P can be sent to the TCM decoder. To be 
able to analyze this receiver, we shall make the simplifying assumption that 
the sequence of states (f>k is delivered to the TCM decoder without errors.^ 

With this receiver structure, we generate (see Fig. 4.1) an equivalent 
discrete-time, discrete- input, continuous-output channel, whose inputs are 
the TCM encoded symbols and whose output is a sequence of quantities 
used to compute the metrics to be associated with the S 3 rmbols This 
channel should resemble a memoryless chstnnel (due to the ideal interleav- 
ing/deinterleaving process) perturbed by additive noise and independent 
fading. Notice that with this transmission scheme the design of the TCM 
scheme is simplified, because we may take adv 2 intage of the wide set of results 
available on the design of TCM schemes for these channels [44]. 

^This assumption is similar to the standard assumption usually accepted in the study 
of coherent detection systems, i.e., that the carrier phase recovery system operates with 
zero error. In both cases this is not realistic, since it implies an infinite signal- to- noise 
ratio, and hence zero error probability. Consequently, the results that follow should be 
interpreted as bounds to the performance of any real-life system. 
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4.2.1 Computation of Rq 

Given a transmission channel to be used for coding, a suitable one-parameter 
characterization of its coding capabilities is provided by the cutoff rate Rq 
created by the modulation system [47]. In fact, although no formal proof is 
available, it is widely believed that Rq is the rate beyond which it becomes 
very expensive to communicate reliably over a channeL 

The computation of Rq is described in the Appendix. For the AWGN 
channel without fading, (4.7) specializes to 


Ro = 2log2iW -log2 


MM 

M + g-|l»(2.-Af-l)-i(2i-Af-l)l|V8A^o 

i=l j-l 


(4.10) 


where in our case 

||s(2i-M-l) - s(2j-Af-l)f 

_ f^‘ II gj2irfc(2i-Af-l)g(t) _ ^j2irh{2j-M-l)q{t) ||2 

Tj Jo 

= ^^T,-2j\os4Trh{i-j)q{t)dt . 

For example, in CPFSK we have 

t 


?(<) = 


2T’ 


and hence 


II s(2i - 1) - s{2j - 1) |p= 4£, 


sin2T/i(t — jy 
2irh{i — j) . ' 


The IRC phase pulse is defined as 


q{t) 


t 1 . „ ^ 

— T*sm2x— , 

2T, 47t T,’ 


and the distances || s(2i - 1) - ■»(2j — 1) || cjin be evaluated by nxomerical 
integration. The corresponding values of Rq are plotted, along with Rq for 
coherent PSK, in Fig. 4.2 for Af = 8 amd in Fig. 4.3 for M = 16. It is 
noticed that for low values of the modulation index h (i.e., small spectrum 
occupancy) PSK has a larger Rq, due to its better distance properties (which 
are traded for the better spectral properties of CPM). 
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When fading is present and perfect state information [64] is available, (4.7) 
becomes 

/ - ^ 


Ro = 2 logj M - log 2 


M M 




V 


j=i j-i 


(4.11) 


Now, by using the results of [64], we obtam, for a Rician channel, 

.2/^. .1 ^ (4.12) 


1 + ^ 

E[«p -f Ail = “P-i + jr + ftJ’ 


while for a Rayleigh channel 

E[exp -p^/?.j] = 

By substituting 

ft, = 8^ II '(2i - 1) - - 1) II’ (■'■I'*) 

in (4*11), we get the desired expression for Rq. The corresponding values of 
Rq are plotted in Fig. 4.4 for M = 8 and Fig. 4.5 for A/ = 16. 

Fig. 4.6 shows the effect over Rq of the selection of the modulation index 
h. If trellis-coded modulation, based on an 7nf{m + 1) convolutional code, 
771 = log 2 Af — 1, is to be used in conjunction with M-ary CPM, we are in- 
terested in the signal-to-noise ratio necessary to achieve a value of Rq equal 
to log 2 M - 1, the actual number of information bits carried by each Ts- 
second signal. This value is plotted versus h for octonary CPFSK. It is seen 
that the lowest values of Ef,/Noj and hence the most power-efficient coding 
channel, is generated for higher values of h* This fact can be explained by 
observing that higher modulation indices correspond to larger phase fluctu- 
ations in the transmitted waveforms, which means that the waveforms are 
more easily distinguishable. On the other hand, larger values of h imply a 
broader power density spectrum, so that the trade-off is between spectrum 
occupancy and power efficiency. This is illustrated in Fig. 4.7, which shows 
the out-of-band power of IRC- CPM as a function of the modulation index. 
It may also be observed that for large h the power efficiency of CPFSK can 
be greater than that of PSK. 

Figure 4.8 shows, for M = 8 and h = 1/4, a comparison between Rq 
associated with optimum and suboptimum coherent detection of CPFSK. It 
is seen that for iZo = 2 suboptimum detection causes a loss in signal-to-noise 
ratio of about 2 dB. 
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4.3 Suboptimum noncoherent detection 

As noted before, to avoid the state growth caused by optimum detection 
and in order to be able to perform interleaving/ deinterleaving, we must get 
rid of the term + <f>o appearing in (4.8). Here we do this by considering 
noncoherent symbol-by-symbol demodulation of CPM signals. 

4.3.1 Computation of Rq 
Let us define the quantities 

A, = (4.15) 

JiT, 

(t+l)T 

A^= / * \{t)s\t,a\)dt. (4,16) 

JiTs 

These quantities, for given values of the fading variable and of the trans- 
mitted symbol aj, are complex-valued Gaussian random variables. Since we 
assuming perfect interleaving, the sequence {A,} is an independent, identi- 
cadly distributed sequence. 

We assume here the following noncoherent detection rule: choose the 
symbol sequence {aj} if 

^|A,f = max^lAip, (4.17) 

i t 

where the TUATimum is taken over all the possible sequences {ui}. 

The pairwise error probability P[a' — + a'] that the detector chooses a' / 
a' when it is known that either a' or a' was transmitted may be upper 
bounded by using the ChemofF boimd: 

P[a' a'|p ] < nunE |exp(-^[^ |A,p - ^ |A,p])|p ,a'| (4.18) 

= ^un^E■[exp(-^[|Aip - |A,p])|pi,o(} (4.19) 

= imn]^E|e^^Mp„7,}, (4-20) 

where ^ is the ChemofF boimd parcuneter (to be optimized to get the tightest 
possible bound) and 

7. = 1A,|2 - lA.f . 


(4.21) 
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The actual computation of the Chernoff bound requires the calculation of 
the expectation E over 7 ,, i.e., over the pair of complex random variables 
Ai, A,. Let Zi denote the following Gaussian random vector: 


Z| — [A|jA,‘], 

so that we can write 

7. = lA,p - |A.l2 = A.A* - A.A* = z*Azf 

where 

A = 


(4.22) 


1 0 
0 -1 


and the superscript T denote transpose. Thus, the computation of the 
ChemofF bound is based on the evaluation of an average of the form 

E |exp ^z" A z^Ip, a' I = J exp ^z*Az^/z(z|p, a')dz (4.23) 

where /z(z|p, o') denotes the conditional Gaussian probability density func- 
tion 


/z(z|p,a') = ^^^^\|Q[ -e3cp{-^[z - z]*Q ^z - z]^} 
with Q the conditioned covariance matrix 


Q=2E 


exp |[z - z]*[z - z]^| Ip, aj 
and z = E[z|p,a']. Explicitly, we must compute the integral 
I - i[z - 
After a cert 2 un amount of algebra, we find 
E|exp[^z*Az^]lp,a'} = 


where 

cind 


B = -iz*Q-‘[I - (I - 2^QA)- V 


IHl 


(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 


IQI |I-2^QA|’ 


(4.29) 
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so that, in conclusion, 


exp{-izrQ-Ml-(l-2eA)-^]}zf 


P[a' ^ a'|p ] < nunn ^ jj _ 2 ^q.a1 


(4.30) 


We now determine the entries of the covariance matrix Q,. Under the as- 
sumption that the noise is white and complex with power spectral density 
2Nq we obtain 

1 A. 


where we define 


Q, =2NoEs 


Ai = -^ [ exp[j2xh{a'i - a')g(t - tT,)]A (4.32) 

T, JiT, 


Ar 1 


and E, is the average signal energy in a T,-second interval. 
We 2 ilso have 

z,' = 2piEs[Ai , 1], 


eind 


|I-2^QiA| = det 


1 - A A. A 
-A*A 1-l-A 


= 1 - 


(4.31) 


(4.33) 


where 

and 


A = 4^NqE, 
Pi = 1 - lA.f . 


After some further lengthy but straightforward computations, we finally 
obtain 


P[a' - s'l p 1 < ■“■‘II '‘‘• 


34) 


It can be easily proved that, under the constraint 


0 < A < 1, 

the exponent of the right-hand side of (4.34) is always negative, as it should 
be to provide a meaningful bound. 
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By using a result in [64], we observe that, for a Rician chainnel, 

2.1 1 + ^ 

l + KT fi "1 + K+/3' 

while for a Rayleigh channel, as A" = 0, 


E[exp -p^/3] = 


Thus, by substituting 

^ No 1- 

in (4.34), we obtain, for any A G (0, 1), 


P[a' 


a'l 


< n 


l^K 


(l-A2/2.)(l + ii:) + fj(A-A>i 


exp< - 


(1 - A2/ii)(l + AT) 4- t^(A — A2)/i, 




for Rician fading, while for Rayleigh fading we find 


P[a' 



1 


(4.35) 


The corresponding Vcdues of Rq are plotted in Fig. 4.9 to 4.12, where A 
WBLS chosen equal to 0.5, a value which was found optimum or near optimum 
in most cases. Fig. 4.13 shows the effect over Rq of the selection of the 
modulation index h. The signal-to-noise ratio necessary to achieve a value 
of Rq equal to log 2 M — 1 is plotted versus h for octonary CPFSK. As for 
coherent detection, the lowest values of JS^/iVo, and hence the most power- 
efficient coding chainnel, correspond to higher values of h. 
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Figtire 4.1: Block diagram of the transmission system. 
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Figure 4.2: Comparison of cutoff rate Rq of coherent octonary PSK, 

IRC-CPM and CPFSK with coherent, symbol-by-symbol detection and ad- 
ditive, white Gaussian noise channel. 
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Figure 4.3: Comparison of cutoff rate Rq of coherent 16-ary PSK, IRC-CPM 
cind CPFSK with coherent, symbol- by- symbol detection and additive, white 
Gaussiein noise chamnel. 
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Eb/No (dB) 


Figure 4.4: Comparison of cutoff rate iZo of octonary coherent PSK, 

IRC-CPM and CPFSK with coherent, symbol-by-symbol detection and 
Rice/Rayleigh fading channel. 
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Figure 4.5: Comparison of cutoff rate Rq of 16-ary coherent PSK, IRC-CPM 
auid CPFSK with coherent, symbol-by-symbol detection and Rice/Rayleigh 
fading channel. 




Eb/No (dB) 



h, Modulation Index 

Figure 4.6: Eb/No for Rq = 2 bits/symbol with octonary CPFSK versus the 
modulation index h. + denote the points at which CPM starts outperform- 
ing PSK. 
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Figure 4.7: Out-of-band power for IRC CPM. The ordinate is the bandwidth 
containing the 99.9% of the signal power. 
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Figure 4.9: Comparison of cutoff rate iZo of IRC-CPM and CPFSK with 
noncoherent symbol-by-symbol detection and of coherent PSK over additive, 
white Gaussian noise chaumel. (M = 8). 
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Figure 4.10; Comparison of cutoff rate Rq of IRC-CPM and CPFSK with 
noncoherent symbol-by-symbol detection and of coherent PSK over additive, 
white Gaussian noise channel. (M = 16). 
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Figure 4.11: Comparison of cutoff rate i?o of IRC-CPM and CPFSK 

w:th noncoherent symbol-by-symbol detection and of coherent PSK over 
Rice/Rayleigh fading channel. (M = 8). 
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Eb/No (dB) 

Figure 4.12: Comparison of cutoff rate Ro of IRC-CPM and CPFSK 

with noncoherent symbol-by-symbol detection and of coherent PSK over 
Rice/Rayleigh fading channel. (M = 16). 
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h, Modulation Index 

Figtire 4.13: Et/No for Rq = 2 bits/symbol with octonary, noncoherent 
CPFSK versus the modulation index h. 





Chapter 5 

Differential CPM 


In this Chapter we describe symbol-by-symbol differential detection of CPM. 
Further, we consider precoding of differentially-detected CPM. The material 
that follows is adapted from [35]. 

The complex envelope of the received signal, say f(t) is delayed by T, 
seconds, transformed into its conjugate f*{t - T,), ^d m^ti^ed by itself. 
Then the real and imaginary components of the signal f(t)f*(t - T,) are 
sampled every T, seconds. As a result, a discrete signal is obtained whose 

phase is 

2X/0T3 + "I” 

where A^„ represents the change over one symbol interval of the signal 
phase, and r;„ represents the change in phase due to the noise. Under the 
assumption that foT, is an integer number, estimate of this phase provides 
a noisy estimate of A^„, which is used to recover the information sequence. 
For fuU-response CPM (the case of precoded partial response will be 

considered later) we have 

Afl„ = Than, 

so that the maximum value taken on by the phase shift is 

= Th(M- 1). 

This quantity must be lower than t because the phases are observed modulo 
2t. Thus, we must choose 
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5.1 Precoding differential CPM 

A form of precoding is developed for use in the channel, to aUow symbol-by- 
symbol detection even for partial-response signalling, in which the frequency 
pulse duration is longer than one symbol time. This precoding undoes, be- 
fore it is done, the correlation among symbols introduced by the inodulator 
and the differential detector. The precoding is based on a modulo-M op- 
eration. Precoding is not needed for fuU-response CPM when the mimm^ 
modulation index is desired; for a larger modulation index or for partial- 

response signalling, it is needed. , i * 

Our description of the precoder will be based on the CPM modulator 
model developed by Rimoldi [34]. The source output (see Fig. 5.1 is a se- 
quence u of independent, equaUy likely, M-ary symbols The first element 

in the transmitter is the precoder, the output of which is the sequence u of 
M-ary symbols u'„. This sequence is input to the continuous-ph^e modula- 
tor. The output of the modulator is an RF signal with information- carry mg 
phase V'(f). The receiver consists of a differential detector and a baseband 
signal processor (BSP). The next element is a phase detector, with output 

AV'(t) + noise, modulo 2x, where AV»(t) = V»(0 ~ ~ 

is sampled every T, to yield the sequence AV»„ + noise, modulo 2X. This is 

input to the BSP, which outputs the sequence of soft-decision estimates of 

the Un. 

We write the CPM signal phase in the form 

exp{j[2T/it-i-V'(t,u')]}. 

In particular, following the notation of [34], we define the “tilted phase” 
by 

V>((r -I- n)T„ u') = 2irh u- -1- irh ^ g((r + + W{rT,) 

i=o »=o 

for 0 < T < 1. Here W is a time function dependent on h, M, and g. The 
frequency /i is lower than the signal’s center frequency by h{M -1)I2T,. 

The differential detector outputs noisy samples of the differential phase. 
Let us suppose that the sampling epoch is r^T,. Without loss of generality, 

T, G (0,1]. 

We define the samples qi of q by 

9t = + 0^*) for t = 0, 1, . . . , L — 1 
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Figure 5.1: Block diagram of differential CPM detection. 


Channel 
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while requiring that 

qo > 0. 

The phase at the n-th sampling time is 
i’n = i’ii'Ts + n)T,yu') 


(5.2) 


= 2 Th’£\' + 4xh^'u;_.g. + Tr(r,T,) for n > 0 

i=0 4=0 

Tlie n-th. sampled, differential phase Ai/^n Is then 
A^„ = V’n - ■0n-I 

= 2x/i [u'„ 2qo + ^ «n-i(23i “ 2«.-i) + «n-t(2 2 ~ 

i=l 

for n > 1. For the case of fiiU- response CPM (i.e. i = 1), we may choose 
T, = 1 so go = 1/2- Then 

AV»n = 2xh u„. (5-3) 

Now we are in a position to see what a precoder can do for us. If 
precoding is not used in the transmitter, so that u' is the channel input 
sequence, then the receiver has to try to remove the effects of Un_u» • • -i^n-L 
from AV-n, the better to detect <. (This is “opening the eye” of the signal.) 
It is better to remove the effects of nearby symbols in the transmitter, where 
the symbols are known without error [32]. To that effect, let us presume 
that u' is the symbol sequence out of a precoder and that u is the M-ary 
symbol sequence input to the precoder and thus to the channel. We want 
AV’n to be closely related to Un only. We have to assume now that all the 
g, ’s are rational. Then 




(5.4) 


do = Q X 2go 

di = Q 'K (2g; — gi-x) for i = 1, . . . , Xi — 1 

di = g X (2 X i -2gL_i) 


where 
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and Q is the smallest positive integer such that all the d^’s are integers. 
In analogy with the coherent reception case [32], we want a precoder that 
makes ^ 

Un = -Rm ( ^ “n-i 
«=0 

where is the modulo-z function. Then we shall have 

Un = RA/(«AV’n), 

where q 

K — _ - • 

2vh 

The precoder has to perform the operation 

L 

v!^ = RjVf («n - ^“n-i ‘^)1* 

i=l 

It is required that do and Af be relatively prime. In the case where L = \ 
and T, = 1, there is no need for a precoder, so u„ = u'„. We can think that 
Q = 1. AV'n and u„ are simply related by 

Un = K. 

The spectrum of the CPM signal is the same with or without precoding, 
since the u'„’s are equally likely and independent when the u„’s axe [33]. 

5.2 Computation of Rq 

In this section we consider the computation of the parameter i?o for the 
channel created by CPM modulation with precoding and differential detec- 
tion. The results obtained from this computation will allow us to choose the 
parameters for the combination of CPM with trellis-coded modiilation. 

We now make two assumptions which allow us to calculate Rq. The 
baseband signail processor (BSP) input (see Fig. 5.1) is R 2 „{Atl>n + noise). 
The first assumption is that the BSP output Vn is simply given by 

Vn = RA/[«R 27 r(At/-„ -|- noise)] e [0,M). 

I 

The second assumption is that 


Q/h = AM 


(5.7) 
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for some positive integer A. The first assumption then takes the form 

Vn = iZM[«(AV»n + noise)]. 

Thus, whether precoding is used or not, the following simple relationship 
between the channel input Un and output V„ can be derived: 

V„ = Rm [Un + «(noise)l (5.8) 

The noise term in the differential detector output, which is the BSP input, 

noise = Pn - Pn-u (^*^) 

where pn is the error in the JV-th noisy phase. It turns out that Rq for 
the channel can be calculated once the probability density function of /i„ is 
known. The pdf of is 


p (a) = J_exp(-7){1 + 2 v^cos(a)[l - Q(v^cosa)] exp(7 cos^ a)} 

(5.10) 

for 0 < a < 2ir, where 

1 r, 

C(a:) = -j= J exp{-t^ I2)dt. 

We can now calculate the cutoff rate Rq for our channel. We derive a 
discrete- input memoryless channel equivalent to our channel in terms of Rq, 
then we present n\imerical residts. From Equations (5.8) and (5.9) we know 
that successive chaumel outputs Vn-i and Vn are given by 


Vn-l = Rm [“n-l + «(Mn-l “ A^n-s)] 

Vn = RMi'^ + ^Pn- Pn-l)]- (5-H) 


This does not represent a memoryless channel. While the chaimel input 
symbols Un— 1 and Un are independent, the noise terms i A*n— 2 and 


fin ~~ Mn-1 not. ^ 

Here we construct a discrete-input memoryless channel, for which Rq 

can be calculated. In constructing the memoryless channel, let us define a 
new sequence of continuous-valued random variables Yn by 


Yn = RM{Vn + -' + Vi + Vo) 6[0,M) forn>l 
The receiver must observe Vq defined by 

Vo = Rm [«(V’o + Mo)]. 


( 5 . 12 ) 
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Given Vo, for each N the set {Yi, ...,Yn} oi random variables can be cal- 
culated from the set {Vi , . . . , V}v} of random variables, and vice versa. We 
can show that the new channel is memoryless. We find that 


n~l 

y„ = (5.13) 

»=1 

from which we can show that the Yn’s are independent. The fact that the 
distribution of Yn depends only on the value of Xn follows from the relation 


Yn = RxiiXn + + I^Pn), 


where the receiver must know (The latter point means that there 

must be phase synchronism at time zero, which can be obtained along with 
symbol synchronism by the transmission of a training sequence before the 

data.) 

Rq is given by 


■ 1 rM / 

Ro = - logj Xf2 / ^ \/piy\k)p{y\i) 

_M Jo ,_Q 


(5.14) 


where p denotes the conditional probability density function for Yi given 
Xi. For a given SNB,, Ro is determined solely by A and M. 

Figure 5.2 shows how Rq varies with A and SNR, for Af = 2 and A = 
1,2,3, 4,8. From this plot (and similar plots that have been obtained for 
M = 3, 4, 8, and 16, we find thet for Rq in the mid-range, having A = 2 
instead of = 1 represents a 4.5-5 dB loss in SNR. S i mila r ly, having A = 3 
instead of A = 2 represents a 3-3.5 dB loss. Figure 5.3 shows how Rq varies 
with M and SNR, for A = 1, Af = 2, 3, 4, 8, 16. 


5.3 Selection of channel parameters 

In this section we will look at how one makes a good selection of the channel 
parameters, which are as follows: 

1. the modulation index h; 

2. the phase pulse q; 

3. the scaled sampling epoch a positive time no greater than 1 (the 
actual szimpling epoch being t,T,). 
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Here we will give examples of the good selection of channel parameters. 
The families of pulse to be treated are named, as is common usage, by the 
frequency pulse and not the phase pulse. The families will be treated m 
order of increasing smoothness, which corresponds to the signal spectnm 
having a wider mainlobe and faster roll-off at frequencies far from the carrier 

frequency. , , . . , 

The highest quedity chaimel (A = 1) can be obtamed with a frequency 

pulse of duration IT if /i = L/M and I is as follows: for the rectangular 
(LREC) pulse, L is any integer; for the triangular (LTRI), L is any mteger; 
for the raised cosine (LRC), X = 1 or 2; for the convolved-raised-cosme 
(LCRC), X = 1 or 2. The second highest quality channel {A - 2) can be 
obtained with h half the size. In general, larger h means wider spectrum, so 
we have Usted only the smallest h in each case. For a given X, different h s 
are provided by different r’s. In general, larger X means narrower spectrum, 
but larger X calls for larger h to maintain the channel quality, so it needs to 
be investigated which X provides the best spectrum. 


Chapter 6 

Doppler phase shift removal 


In this Chapter we consider the effect of a Doppler frequency shift over 
differentiaUy-detected CPM signals, and we derive a circuit that makes it 
possible to estimate and to remove the corresponding phase shift. 

Let us write the observed signal r(t) in the form 

( 6 . 1 ) 

where P{t) and v{i) are the amplitude fluctuation and the phase fluctuation, 
respectively. They account for fading, noise, and intersymbol interference, 
/d is the Doppler frequency shift. 

For differential processing, the signal (6.1) is synchronously sampled at 
times 

in — n = 0) 1) • • • f 

and the following sequence is formed: 

( 6 . 2 ) 

A 


where 

4>d = 

is the Doppler phase shift, 

A6n = ^(tn,a) - d(tn-i,a) = vhan 
represents the phase fluctuations due to the transmitted data, and 

r]n = »'(in) - t'itn-l)- 
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FinaUy, 

Pn = 

(Perfect symbol syncbronization is assiimed.) 

It is seen from (6.2) that the presence of a Doppler frequency shift causes 
the phase of the sequence fn to be altered by a term <j>d, added to the 
information sequence This has to be removed prior to demodulation. 


6.1 Estimation schemes 


Two constrzdnts are associated with the design of a Doppler compensation 
circuit for continuous-phase modulated signals with application to mobUe 
satellite communications: 


• F 2 ist frequency ^lcquisition is required. In fact, if data are transmit- 
ted in short bursts or packets, the acquisition time should not be a 
significant portion of the burst interval. 

• Since multipath fading affects the propagation, the Doppler estimator 
must be insensitive to the signal amplitude fluctuations caused by 
fading. 


As a consequence of the first requirement, we consider open-loop esti- 
mation structures, as first suggested (for different systems) in [82,78,29]. 

The class of methods we propose to estimate the Doppler phase shift 
<f>d is based on the following considerations. Consider a random VMiable 0, 
observed mod 27 t/A and taking values 6 e (-x/A,)r/A). Let f{e^ ) be its 
probability density function. We assume that /(•) is a symmetric function. 


SO that 


J-2ir/A 


is a real quamtity. Now the problem is the following. Assume that we observe 
K (say) values of $ = 0 + where is a constant, and we want to estimate 
<f> based on these observations. The probability density of 0 -b ^ is obtained 
by shifting (mod 2ir/A) the original function f{6), so that we get 


J-2ir/A 
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Thus, under our symmetry assumption, <f> can be obtained from the phase 
of E[e-'^®], which in turn can be estimated as 

] Si ^ (6-3) 

" fe=i 

where V’fc, * = 1, . . . , if, denote the observed values of 

Observation of eq. (6.2) shows that the phase of the received samples 
includes a noise term which we shall assume to have a symmetric prob- 
ability density function, and a data term A0„, which takes on values 

±rh,±37Th,...,±(M -l)vh (mod 2x). 

Since h = JJM, the probability density function of the above phase turns 
out to be periodic with period 2wh, and in the interval (— t,t) it has copies 

centered at iirh, i = ±1, ±3, . . . , ±(M — 1). 

From the theory presented before we see that the estimate of the Doppler 

shift (f>d can be based on the phase of the quantity 


zn = ( 6 - 4 ) 

and the positive function F(-) is arbitrary. Thus, the estimator of the 
Doppler phiise shift becomes, if is a positive quantity, 

^ " M j. eS F(j>n) co» + A«„ + .h) 

In words, (6.4) is equivalent to performing a transformation from rectangular 
to polar coordinates on each complex sample r„. Next, we perform two 
transformations on its amplitude and phase. Finally, we perform a polar-to- 
rectangular transformation on the result. As observed in [82], in a practical 
implementation the nonlinearity becomes a read-only memory transforming 
a quantized complex number into another quantized complex number. 

The choice of the ftmction F should take into account the presence of 
fading. In fact, samples heavily affected by fading should be given a lower 
weight in the time averages which form the estimate In our simulations, 
we have always assumed F{pn) = Pn, which is perhaps the simplest choice to 
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decrease the influence of the samples whose amplitude is lowered by fading 
cflfect s« 

For future reference, observe that (6.5) can be rewritten in the form 

( 6 . 6 ) 


1 Ac sin M<j>d + -^5 cos M<f>d 

As cos M<f>d — Ac sin M<l>d' 


where 


K-l 


Ac = cos Af(A0n + ^n) 

n=0 

= 4 X! sin Af(A^„ + »7n) 

n 


n=0 


The trigonometric identity 


, jB sin V> + A cos il> -i A 

tem ^ -= ; — 7 — ^ + tan , 

5cosV'-AsmV’ 


(6.7) 


in conjunction with (6.6), provides an expression for the estimation error: 

( 6 . 8 ) 


1 j. L . -i As 


0.1.1 Estimator A 

The first estimator we consider is based on eq.(6.5), with F{pn) = Pn- 
6.1.2 Estimator B 

We consider here an estimator structure whose aim is to broaden the range 
of Doppler phases that can be tracked. 

Assume that the received signal is sampled twice per signaling interval, 
namely, at time instants (n + a)T, and (n + b)Ta, n = 0, 1 , . . ., where 0 < 
a < 6 < 1. If one sample is multiplied by the conjugate of the other one, we 
get a sequence of quantities whose phases have the form 

= (6 - a)<f>d + 2Tha„A(a,6) + q„, 

where A(a, b) < \ denotes the area of the frequency pulse used for CPM 
between times aT, and bT,. (Note that q„ has now a meaning different from 
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before). The presence of the factor 6 - a in the Doppler phase term causes 
an increase in the estimator variance: in fact, we have 

instead of (6.8). Under the rough assumption that the term tan"! ^ re- 
mains the same, the estimation error increases by a factor l/(b - a). This, 
in turn, increases by a factor of approximately 1/(6 - af the value of the 
sample size K necessary to achieve a given estimator variance. 

As a consequence of this estimator structure, the probability density 
function of the phase turns out to be periodic with period 47rM(a,6), a 
quantity less than 2^h. To avoid ambiguities in the estimation procedure, 
we must have: |(6 - a)M < irhA(a,b). Thus, with respect to Estimator A, 
we get a broader range of Doppler phases that can be estimated if the ratio 
p = 2A{a,b)lb - a is greater than 1. 

0.1.3 Estimator C 

This estimator further extends the range of Doppler frequencies that can be 
1)T£tclC€d« 

The basic idea underlying Estinoator C is that in one symbol mterval 
the phase variation due to the Doppler frequency shift is linear, while that 
due to data modulation depends on the shape of the phase pulse, and hence 
can be made nonlinear. Thus, it exploits the nonlinearity of the phase pulse 
by performing two differential detections based on four samples within each 
symbol interv 2 d. Details about it can be found in [36]. 

6.2 Simulation results 

We now describe some results arising from computer simulation of Doppler 
frequency shift estimators A,B, and C. We assume an octonary, full-response 
CPM with modulation index h = 1/8 and IRC frequency pulse, i.e., 

1 , 2irt . 

In Estimator B, a = 0.1 and h = 0.9, so that 
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and 


P = 


2A 

b - a 


= 1.234. 


The comparisons are made in terms of the variance of the Doppler frequency 
estimation error. For all the estimators, 4095 symbol intervals were simu- 
lated. For Estimators A and B, a window of K intervals was used to compute 
(^, and hence /j. This window was slided by one symbol interval at the time 
to generate a set of estimates. The estimate of was then obtained by av- 
eraging the squared estimation error over 4095 - K runs. (Notice that with 
this procedure the estimates become less and less accurate as K increases). 
For Estimator C, the time-average was evaluated by using a one-pole But- 
terworth filter whose output was sampled at time KTj. 

In Figs. 6.1 to 6.6 the vari^mce of the Doppler frequency estimation error 
is plotted vs. K for several vedues of signal-to-noise ratio. Figs. 6.7 to 6.9 
compare the error variances of estimators A,B, and C for three values of 
signal- to- noise ratio. 


6.3 Conclusions 

Three estimators of Doppler frequency shift in CPM signals were proposed, 
and their performance evaluated by simulation. Their structure is based on 
differential detection, which makes them attractive for use in a differentiaBy 
coherent receiver (although in principle they can be used in conjunction with 
any other CPM detection scheme). 

These estimators are in order of increasing complexity. Increased com- 
plexity is made necessary to remove ambiguities in the estimate when the 
Doppler frequency shift is high. Moreover, at least for high signal-to-noise 
ratios, a wider range of tradeable Doppler frequencies causes a substantial 
increzise in the variance of the estimation error. 



ERROR VARIANCE 


CHAPTER 6. DOPPLER PHASE SHIFT REMOVAL 


PERFORMANCE OF DOPPLER ESTIMATOR 3 

10000 
1000 
100 
10 
1 

.1 

.01 
.001 
.0001 

0 2 4 6 8 1 0 12 1 4 1 6 1 8 20 22 24 26 28 30 32 34 36 38 40 



Figure 6.1: Estimator A. Variance of the Doppler frequency estimation error 
in the presence of AWGN and data modulation, fd = ^ Hz, vs. the number 
K of samples averaged (abscissa label is iiT/lOO). 
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Figure 6.2: Estimator A. Variance of the Doppler frequency estimation error 
in the presence of AW^GN and data modulation, = 134 Hz, vs. the number 
K of samples averaged (abscissa label is FT/lOO). 
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Figure 6.3: Estimator B. Variance of the Doppler frequency estimation error 
in the presence of AWGN and data modulation, /(i = 0 Hz, vs. the nTimber 
K of samples averaged (abscissa label is iiT/lOO). 
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PERFORMANCE OF DOPPLER ESTIMATOR 4 



Figure 6.4: Estimator B. Variance of the Doppler frequency estimation error 
in the presence of AWGN cind data modulation, = 134 Hz, vs. the number 
K of samples averaged (abscissa label is K/IQO). 
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Figure 6.5: Estimator C. Variance of the Doppler frequency estimation error 
in the presence of AWGN and data modulation, fd = 0 Hz, vs. the number 
K of samples averaged (abscissa label is K/lOO). 
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PERFORMANCE OF DOPPLER ESTIMATOR 6 



Figure 6.6: Estimator C. Variance of the Doppler frequency estimation error 
in the presence of AWGN and data modulation, fj = 134 Hz, vs. the number 
K of samples averaged (abscissa label is iif/100). 
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COMPARISON OF DOPPLER ESTIMATOR 
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Figure 6.7: Variance of the Doppler frequency estimation error in the pres- 
ence of AWGN auid data modulation, SNB,=10 dB, = 134 Hz, vs. the 
number X of samples averaged (abscissa label is if/lOO). 
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Fieuie 6 8: Variance of the Doppler frequency estimation error in the pres- 
ence of AWGN and data modulation, SNR=20 dB, /a = 134 Hz, vs. the 
number K of samples averaged (abscissa label is iif/lOO). 
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Figure 6.9: Varizince of the Doppler frequency estimation error in the pres- 
ence of AWGN and data modulation, SNR=30 dB, fd = 134 Hz, vs. the 
number K of samples averaged (abscissa label is RT/100). 



Chapter 7 


Timing recovery circuit 


This Chapter presents analysis and simulation of a circtiit designed to re- 
cover the symbol timing from a CPM signal. This circuit is based on the 
inherent self- synchronizing capabilities of CPM, and in particidar on the 
property that, when raised to a suitable power, a CPM signal with rational 
modulation index wiU exhibit a power spectral density with spectral lines 
spaced by the clock frequency. 

After the relevant mathematical analysis, we describe the simulation 
procedure. Simulation results are presented and interpreted, with the aim 
of obtaining design guidelines. The conclusion is that, with a proper choice 
of the filters and of the CPM pulse shapes, this circuit can actually solve the 
problem of recovering the symbol timing in our system with a high enough 
accuracy. 

7.1 Generalities on symbol synchronization 

In this section we show the basic operation of the synchronizer. Later on, 
we shall expand on this analysis to illustrate how to modify the synchronizer 
to allow it to work in the presence of fading and Doppler frequency shift. 
The original design was suggested in [3], which presented the system as a 
general method for recovering symbol timing and carrier wave phase and 
frequency in a CPM system- For our purposes, the system was simplified to 
recover only the symbol timing. Fig. 7.1 shows the elements of the system, 
with typical filter bandwidths for an 8-ary CPM system, when = 2400, 
h = 1/8, and in the absence of Doppler shift. 

The analysis in this section is the basis for a computer simulation of the 
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system. The graphs in this chapter were obtained from simulation, using the 
above mentioned pairameters and a pulse shape of 1-REC. Later, simulation 
will be used to cedculate the variance of the s3rmbol timing jitter in the 
presence of Rician fading, Doppler shift, and Gaussian noise. 

The motivation behind the system stems from the fact that a CPM signal 
with an integer valued modulation index h contains spectral lines [2], which 
can be extracted to produce the symbol timing information. When Doppler 
shift is present, the lines shift in frequency, but, since they all shift the same 
amount, the timing information is not lost. This will become more clear as 
the analysis progresses. The underlying principle of the system is based on 
two standard trigonometric identities: 

cos^ * = i + ^ cos 2® (7-1) 


cos X cos y=\ cos (* + y) + r cos {z — y) 

2 » A 


(7.2) 


These two identities suggest that, given a CPM signal with modulation 
index, h — I/7, we can change it into a CPM signal with h = 1 by raising 
the signal to the 7'^^ power. This is Ulustrated below, with a simple cosine 
function with argument 


cos^ 0 i + ^cos2fi (7.3) 

3 1 

cos^ d = — cos 5 4- T cos 35 (7.4) 

4 4 

3 1 1 

cos'* 5 = - + - cos 2tf + - cos 4t> (7.5) 

8 2 8 

So, if we pass a CPM signal with h = I/7 through a 7-power device, the value 
of ^ in the above equations is 2irfot + <f>{t,a). The output from the device 
contains several cosine terms, including one with argument 'y2irfot+'Y<f>{t, a). 
Now, since the originsd a) is calculated using h as the modulation index, 
the “new” a) is essentially the same as the original one, except that the 
“new” h is equzd to 7 times the original h, or hnew — Unoriginal — 1* In [2] 
it is shown that 2in integer-/i CPM system contains spectral lines, and the 
“new” CPM signal, with h = 1, is in that category. 


7.1.1 Basic Synchronizer 

Figure 7.1 is the block diagram of the basic synchronizer. The operation 
is straightforward. First, BPFl p£isses the received signal, while limiting 
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ti. mount of noise. Next, the CPM signal is passed thiough a j-power 
device, which changes it to integer-/. CPM. Now the onl, 
einansion that is useful is the highest frequency term (q/o), »» BPF2 selects 
tS rone of frequencies near y/o. Next fUter, NBFl NBPyaa 

a spectral line, one located at q/o + n«./2 “<* ~ ^ 

appear in the newly formed integer-/. CPM. The two selected hues m^^en 
multiplied together to form sum and difference ftequenaes The Terence 
frequency is then selected by IBPP, yielding a smusotd of frequency nff.^ 
The frequency of the recovered sinusoid is then div.ded hy n, yieldmg 
signal with frequency the symbol rate. 


7.1.2 Low-Pass Equivalent Signals 


In order to simulate the system described in Fig. 7.1, low-pass eqmvalent 
signals representing the complex envelopes of the bypass «g-als were used 
tlnoughout the simulation. (This is also a possible means of actual hard^ 
ware Splementation.) The advantage of using low-pass eqmvalents mstead 
of their bandpass counterparts is that a lower sampling ra e co e us , 
without any loss of information. This keeps simulation run tunes to a mm- 
imum. For this reason, the frequency spectrum plots shown m 
are for the complex envelopes of each signal. The spectrmn of the actual 
bandpass signals are derived by shifting the spectrum of the coinplex enve- 
lope up in frequency, such that the dc (zero-frequency) component is located 
at the reference frequency of each plot, and each plot, bemg symmetrical, is 

reflected across its own reference frequency. . , , 

Although the analysis below is general enough to mclude the effects 
Doppler shift, noise, and fading, the spectra illustrated in tWs chapter were 
calculated without adding any of these parameters m order to Jl^^trate 
the basic principles of the system without any complications The effects 
those parameters on symbol timing recovery will be mcluded later. 

7.2 System Analysis 

7.2.1 Signal Definitions 

We begin the analysis by referring to the definition of the CPM sign^. 
Th. spvctrmn of this sigoal (point “A” on Pig. 7.1) is shown m Fig. 7.2 
This piot is tho spectmm of the complex envelope of the s.^M, which, 
teferemed to the cnrriet frequency of the CPM signnl. This spectmm 
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was calculated using the computer simulation to provide a time-domain 
representation of the signal, and then taking the discrete Fourier transform 
of the data points. The sampling frequency of the simulation was /, = 16iZj. 
BPFl was designed to have a bandwidth sufficiently wide to pass the CPM 
signal undistorted, and to allow for a given amount of Doppler shift, while 
limiting noise. Therefore, the spectrum at point “B” is identical to the 
spectrum at point “A”, except that the noise power at “B” is reduced. A 
theoretical power spectrum is shown in Fig. 7.3, which was calculated using 
a technique described in [9] cind presented in greater detail in what follows. 
In Fig. 7.3, the frequency rainge is normalized to the symbol rate, iZ,, but 
the plot is otherwise in agreement with the simulated results. Notice that 
the magnitudes of the two power spectra are different, but that the relative 
magnitudes are essentially the saime. There are two reasons for the difference 
in absolute magnitude. First, the theoretical graph is normalized to the 
symbol rate, so the magnitude is scaled to the normalization factor. Second, 
the simulation results were calculated using a discrete Fourier transform, 
so the resulting power level depends on the block-length N. Notice that 
the TTiaiu lobe of the signal is fairly narrow, with a low-pass-equivalent 3 
dB bandwidth of about lZ,/4 (the actual band-pass signal has twice this 
bandwidth, or Ra/2). 

At the output of BPFl, we have the undistorted CPM signal with 
Doppler shift. A/, plus bsmdpass Gaussian noise, n(t) (In this section A/ 
and n(t) are both set to 0, but they are included for completeness.) We let 


the signal power be P. Thus 

y{t) = cos{2v{fo + ^f)t + + n{t) (7.6) 

= I{t) cos 2ir fot — Q(t) sin 2irfot (7-7) 

where 

I{t) = cos (2xA/t -f a)) -f Tic(t) (7.8) 

Q{t) = sm{2ir A ft -f ^(t, a)) -|- n,(t) (7.9) 

and Tic(t) and nj(t) are the in-phase and quadrature-phase noise compo- 
nents, respectively. In polcir form: 

y{t) = R{t) cos (2T/ot -f- 0{t)) (7-10) 


R{t) = ^P{t) + Q^t) 


0(t) = tan ^ 


m 

m 


where 


(7.11) 

(7.12) 





CHAPTER 7. TIMING RECOVERY CIRCUIT 


105 



7.2.2 7 -Power Device Operation 

At the output of the 7-power device, the signal enters BPF2, which serves 
as a zonaJ filter, allowing components in the vicinity of 7/0 to pass through 
it. The signal in the vicinity of 7/0 is: 

y^(t) = [il'^(t)cos^(2x/ot-|-0(t))].,/o 

= cos (7^ (t)) cos 2x7 /ot — sin (7^(t)) sin 2x7 /ot] 

= /-y(t) cos 2x7 /ot - Q-y(t) sin 2x7 /ot 

= Re[iT(t)e^'^'’(‘)e-''2’^^/o‘] (7.13) 

where ib is a constant, and equation (7.13) is a bandpass signal with complex 
envelope (referenced to 7/0): 

Uit) = (7.14) 

Without noise, and with h = I/7, 

j 00 00 

7^(0 = = 7(2x-) ^ aiq(t - iT,) = 2x ^ a,g(t - iT,) (7.15). 

7 i=_oo i=— oo 

which corresponds to integer-valued h CPM, having spectreJ lines located at 
odd integral multiples of Rs/2^ as shown in [2]. The spectrum at point 
is shown in Fig. 7.4, which is the low-pass equivalent spectrum referenced to 
fo- The actual bandpass signal has 8 spectral lines centered about /q. The 
spectral lines are very striking in that plot, confirming the results in [2]. 

7.2.3 Filters NBFl and NBF2 

Next, we need to select certain spectral lines using bandpass filters NBFl 
and NBF2. The choice is somewhat arbitrary. In fact, any narrow filter 
should work (without Doppler shift). Since the signal is defined as a low- 
pass signal, the bandpass filters must also be described in terms of their 
low-pass equivalents, with in-phase and quadrature-phaise impulse response 
components of each filter. 

We will analyze the case for NBFl; the analysis for NBF2 is identical, 
except for several obvious substitutions. Starting with the signal entering 
the filter, we realize that it is a low-pass signal referenced to 7/0. NBFl 
has a center frequency of (7/0 + nRs/2)^ which means it caniDe described 
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by a low-pass filter, referenced to its center frequency. It has an impulse 
response: 

h{t) = (7.16) 

h{t) = h,{t) + jh,{t) (7.17) 

and hc{t) emd ha{t) are the in-phase and quadrature-phase impulse responses 
of the filter. Rewriting the bandpass signal, x{t), entering NBFl as: 

x{t) = Re[i(t)e^'2"^'f“‘] (7.18) 

= Re[i(t)e"^'^"?^'*e-'^’"<^*'^^^*^‘] (7.19) 

= Xc{t) cos 2 t( 7 /o -I- ^Rs)t - x,(0 sin 2 t( 7 /o -I- (7.20) 

wli6r€ 

i{t) = Xc{t) -1- jx,(t) (7.21) 

Xc{t) = Re[i(t)e--'2’^T«-‘] 

= Xc{t) COS (27r^R,f) -i- *,(t) sin (2ir^R,t) 

A « 

x,(t) = Im[i(t)e-^2’^5^*‘] 

= -Xc{t) sin (2ir^R,t) -|- *,(t) cos (2x^R,t) 

where Xc{t) and *,(0 are the in-phase and quadrature components of the 
low-pass equivalent of x{t). 

Now, with this new signal referenced to the same center frequency as 
NBFl, we use the equivalent low-pass filtering technique described in [10] 
and diagramed in Fig. 7.5. The digital filters used were designed using At- 
lanta Signal Processors Digital Filter Design Package Version 1.02. With 
symmetrical bandpass filters referenced to their center frequencies, only the 
in-phase component, hc(t), of the impulse response is needed, since the 
quadrature- phskse component, is exactly zero. Therefore, Fig. 7.5 

simplifies to Fig. 7.6. A similar analysis holds for NBF2, with x(t) being 
re-referenced to (7/0 — nRs/2). Thus, we produce yi(t), the filtered version 
of x{t) through NBFl, and ^he filtered version of r(t) through NBF2. 
The power spectrum at point “D” is shown in Fig. 7.7, and at point “E” 
in Fig. 7.8 using n = 1. In Fig. 7.7. These two figures are plotted slightly 
differently from the others. They include both the positive and the negative 
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frequencies of the spectra, since the spectra are not synunetrical about their 
center frequencies, and it is importaint to show what is occurring on each side 
of the center frequency. In Fig. 7.7, the center frequency is (7/0 + 1/2T,). 
Notice that there are 3 additional spectral lines that are higher in frequency 
than the selected line, and 4 lines lower in frequency. Although the presence 
of these additional spectral lines is visually apparent in the graphs, they are 
attenuated about 55 dB from the spectral line of interest, due to the atten- 
uation of the narrow-band filters. Since they are greatly attenuated, they 
actually have little effect in the spectrum. In Fig. 7.8, the center frequency 
is _ 1/2T,). In that graph, there are 4 spectral lines in the positive 
frequency direction, and 3 lines in the negative direction. The difference 
between the two plots is related to the spectral line that we are recovering; 
the first one has 3 lines above it, and the second one has 4. But in either 
case, the additional lines are greatly attenuated, so as to be insignific^t. 

In Figs. 7.7 and 7.8, there are basically two discrete frequency signals 
(one at 7/0 + l/2T„ and one at 7/0 - 1/2T,) plus some continuum, due 
to the randomness of the data stream and any additive noise in the signal. 
(If the data stream were periodic and there were no additive noise, there 
would not be ainy continuum.) At point “D”, the signal is a sinusoid of 
frequency (7/0 -l-nff,/2). At point “E”, it has a frequency of (7/0 - nRj2). 
When the two are multiplied together, they produce a signal that has the 
sum and difference of the two frequencies. The resultant waveform contains 
three types of components: line x line (the desired components), line x 
continuum, and continuum x continuum. Obviously, the first component is 
what we need to produce our final results, while the effects of the additional 
components must be minimized. 

7.2.4 Multiplier and LBPF 

Next, we multiply these two bandpass signals together, and keep them in 
low-pass terms, referenced to a common center frequency. We have the 
two signals, yi(t), which is referenced to (7/0 -I- niZ,/2), and j/2(t)> which is 
referenced to (7/0 - riR,/2). Writing them as bandpass signals, we have: 

yi(t) = Vci (0 cos 25t(7/o 4- ^-K»)t - y*i (*) 2t(7/o -I- (7.22) 

y2(t) = ycj(t) cos27t(7/o - ^Ra)t - Va:iit)sm2'ir{'ffo - (7.23) 
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Multiplying them together yields 

yi(0y2(0 = [cos 2x27 /ot + cos 27rnjR,t] (7.24) 

— (0 2Tr2^fot — sin27rniJ,t] 

2 

_yc2(0yii(0rgjj^2xniZ,t + sin 2x27 /ot] 

2 ^ 

y5i(0y^2 (0[cos2yn.R,f - cos2x27/ot] 

2 

which contains frequency components around 27/0 and niZ,. This suggests 
that we can now recover either the carrier frequency (by high-pass filtering 
and dividing by 27), or the symbol timing information (clock rate). To 
recover the symbol timing, we select the components around nRf. In the 
vicinity of niZ,: 

^(^) = yi(0y2(0 = ^(yci(*)yc2(*) + y*i(0y»2(*))cos2xniZjt (7.25) 
-\ivc2 (<)y.i (0 - y^2 («)yci (0) si® 2xniz,t 
So, the low-pass eqmvalent signal referenced to niZ, is: 



z{t) = Re[i(t)e^'2’^"^*‘] 

(7.26) 

with 


z{t) = Zc{t)+jZs{t) 

(7.27) 

and . 


2c{t) = -(yci(0ycj(0 + y*i(0y*j(0) 

(7.28) 


^s(0 ~ 2^^‘^2(0y*i(0 ~ yij(^)yci(0) 

(7.29) 


The spectrum of this recovered signal is shown in Fig. 7.9, representing the 
signal at point “F”. In the simulation we have let n = 1, so the spectriim in 
the figure is referenced to nR, = iZ,, the symbol rate. This signal is then 
filtered by LBPF in a similar way as in NBFl and NBF2. Now, we have a 
signal T(t), which is the filtered version of z{t), with an in-phase component, 
Tc(t), and a quadratme-phase component Tj(t). T{t) will become the symbol 
timing signal, Taymit). The bandpass signal after LBPF is now: 

T{t) = Re[(T^(t) 


(7.30) 
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where 

T 

<f>{t) = arctan(y^j^j) 

(7.31) 

7.2.5 

Symbol Timing Recovery 



In order to recover the symbol timing signal, the hardware would count zero 
crossings of this bandpass signal, and then divide the frequency by n. For 
simulation purposes, we form the signal 


T,ym{t) = cos + <^(0] (7-32) 

which is a sinusoid of frequency 2irRsy and phase, 4^{t)lny that varies slightly 
with time, otherwise known as phase jitter. Ideally, <f>{t)ln would be con- 
stant at a value proportional to the time delay between when the signal 
is transmitted and when it is received. This allows the system to track 
the CPM signal, given varying delay times. Of course, in practice, the 
phase is not constant, due to the noise, fading, and signal dependent jitter. 
Therefore, a suitable performance measure is the variance of <f>{t)y which 
is a measure of how varies with time. The spectrum of the recovered 
signal is shown in Fig. 7.10, representing the signal at point ”G”. Since 
this is the signal of interest, we have chosen to plot the spectrum of the 
actual bandpass signal, rather than the spectrum of the complex envelope 
of the signal. Therefore, the line at Rs = 2400 contains the symbol timing 
information (frequency and phase). The phatse of this recovered signal is 
shown in Fig. 7.11. The variance of the recovered phase in the figure is 
1.3 X 10“^, which is quite small. A more detailed discussion on interpreting 
the phase- jit ter variance is included later. 

7.2.6 Fading Description 

In order to simulate the fading effects on the CPM signal, a Rician fad- 
ing model was used. The complex fading envelope has an in-phase and 
quadrature-phaise component given by: 

^,(0 = (7.33) 

F,{t) = ^s(t), (7.34) 

where A is the direct component, and «;c(t) and «J,(t) combine to give the 
scattered component of the signal. The ratio Pdirect! Pacattered is called K, 
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where 


K = 


A ^/2 


scat 


(7.35) 


and cr^cat power in the scattered component. The scattered component 

of the signal is calculated by first generating two independent sequences of 
white Gaussian variables, each with a variance <rf. The two sequences are 
then filtered by a low-pass filter, yielding a variance (coming out of the filter) 
of which is equal to the power in the scattered component of the noise. 
The bandwidth of this filter is chosen to equal the Doppler spread^ Bd, of the 
channel [17]. The reciprocal of Bd is approximately equal to the coherence 
time of the channel, which is a measure of how fast the envelope of the fading 
signal chzmges. A large coherence time means that the fading signal is slow 
to change, and that the bandwidth, Bd, is fairly narrow. In our simulation, 
we used a slow-fading approximation, which means Bd ■<C Rs (in onr case, 
Bd ^ 100 Hz), and the fading variable was updated only once every symbol 
time, so it stayed constant from symbol to symbol. The relationship between 
the input vari£ince zuid output variaince of a Gaussian random process going 
through a linear, shift-invariant filter is: 



where f, is 
bamdwidth 
given by: 


the sampling frequency, and B is the one-sided noise-equivcdent 
of the low-pass filter (100 Hz in our case). Therefore, is 


<T 


2 _ 
5cat 



(7.37) 


Combining (7.35) and (7.37), we find that <r? is given by 



A^fs 

\KB 


(7.38) 


Once the fading envelope is generated, the original CPM signal is multiplied 
by the fading envelope to produce: 



= Re{(7(t) -b jQ{t)){Fc{t) + jF 3 (t))} 

(7.39) 

Q/atie(0 

= lm{{I{t) + jQmFc{t) + jF,{t))} 
= I{t)F,{t)-{-Q{t)F,{t) 

(7.40) 
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7.2.7 Noise Description 

Finally, in order to add bandpass noise to the signal, we use a low-pass 
eqmvalent noise signal. We must first determine a relationship between <rf , 
the input variance of the Gaussian noise source, and Et,/No, the bit energy 
per noise power density ratio. For a bandpass filter with bandwidth 2B, we 
know the following: 

S P Ei)R^ Rb , Eb \ /y ^1 ■! 

N ~ ~ Nq 2B ~ 2B^ No’ 


Therefore, 


N = 


P2B 

RbiEblNo) 



Setting equation (7.36) equal to (7.42) gives: 




Pfs 


Rb{Eb!No) 


(7.42) 


(7.43) 


Therefore, in order to simulate a given EbfNo, we generate two independent 
sequences of white Gaussian variables with variance <rf, defined by equation 
(7.43). Then, we filter those sequences through BPFl, using the technique 
in Fig. 7.6, 2 md the resulting <7 q produces the desired Eb/No, when the noise 
is added to the signal. A plot of the noise spectrum is shown in Fig. 7.12, 
for EbINo = 10 dB, with a low-pass equivalent bandwidth of 2500 Hz for 
BPFl. The sidelobes in the noise spectrum are a result of using an elliptic 
filter design for BPFl, since the noise spectrum is essentially the frequency 
response of BPFl. 

The nei Section will address the additional problems encountered when 
the CPM signal gets faded, noise is added, and Doppler frequency shift is 
introduced. 
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7.3 Application to Mobile Environment 

The basic synchronizer system described in the last section was designed for 
a non-mobile (stationary), non-fading environment. However, in a mo e 
environment there are additional problems that must be taken mto account, 
in order to ensure that the synchronizer will work. 

7.3.1 Time Delay 

The basic problem that must be solved in any environment is to recover 
the phase of which is proportional to the time delay between the pomt 
when the CPM signal gets transmitted and when it is received, 
mobde environment, accurate phase recovery is import^t ^ 

unknown distance (hence unknown delay time) between the ve^c e and the 
satellite. As an example, a time delay that is ^ mteger m^tiple of T,/4 
would yield a recovered phase, <j>, equal to 2ir/4 - ir/ . ecause o 
extreme importance in recovering this exact phase, the synchronuer must 
be able to unambiguously track the symbol timing over the range 0 to Ts, and 

Notice that the last block of Fig. 7.1 is a divide-by-n function. The 
analysis in the last section suggests the reason why the 

BasicaUy, the frequency of the recovered signal entenng the divide-by- 
block is n times the symbol rate (nR,). Because the p W of the s^ 
gets divided by n as weU as the frequency, a value of n that is anythmg 
other than 1 will produce an n-fold ambiguity in the recovered phase of the 
signal. Take, for example, the case when n = 3. Assume the true P^ase o 
the CPM signal is 2it/2. When multipUed by 3, it becomes 2 t, and after 
the divide-by-n, it becomes 27 t/ 3. Next assume the true ph^e is 4x/3. 
When multipHed by 3, it becomes 4ir, which is indistingmshable from 2x. 
When passed through the divide-by-n, it becomes 2x/3, w^ch is obviously 
incorrect. Similarly, a CPM signal with a true phase of 2x becomes 2x/3 at 
the output of the synchronizer. Clearly, this produces an ^biguity m the 
recovered phase. The result of this ambiguity is the mabihty to track the 

recovered phase for varying signal delays. 

There are several approaches to resolving this ambigmty problem. One 
approach is to recover two pairs of spectral lines, such as the n — 5 an t e 
n = 7 lines. That will produce two phases, <^s(t) and <^7(0» which may e 
different from one ainother: 

Mi) = Mi) + «5 Y 


(7.44) 
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Figure 7.8: CPM Spectrum (Point “E”) 
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Mt) = Mt) + cry . 

The inherent ambiguity may be resolved by finding the smallest mteger 
values of 65 and er in the range (0 < ej < 5; 0 < C7 < 7) such that 

Mi) ^ Mi) 

Aside from the fact that there would still be the possibihty of an e^or 
in the recovered phase, this technique is inherently complicated and time 
consuming. A better technique, in the sense that it is simpler and less prone 
to errors, is to avoid the ambiguity issue completely, by choosmg n - 1. 

The selection of n = 1 can cause difficulty in the recovery process.^ In [4], 
it is stated that the value of n which produces the least phase jitter is 
n = M-1, where M is the alphabet size. If M > 2, then the best value 
of n is greater than 1. However, in resolving the ambiguity issue, we ha,ve 
already chosen n = 1. Therefore, we must be careful, when we select a pu^e 
shape, to choose one that wiU produce as small a phase jitter as possible, 
given that we are using the n = 1 spectral line. This means we need a pu^e 
shape with a strong spectral line at n = 1. Next section includes the results 
that support our preference for 1-RC. 


7.3.2 7 -power device effects 

The system considered contains a component that magnifies some of the 
disturbances (noise, fading, Doppler) that occur in a mobile-satellite en- 
vironment. The component is the 7-power device, which raises the CPM 
signal to the 7 power. Breaking the action into magnitude and phase re^- 
tions, it raises the magnitude to the 7 power, and multipUes the phase by 
7. Therefore, disturbances that affect the ampUtude of the CPM signal are 
exponentially-magnified, whUe distortions that affect the phase of the CPM 

signal eire lineairly-magnified. . .1, 

The negative effects of the 7-power device are especially noticeable on tne 
Doppler shift and the fading of the CPM signal. First, the device multiphes 
the received Doppler shift by 7, making Doppler uncertainty 7 times larger 
than originally anticipated. The two filters (NBFl and NBF2) must be 
chosen to have adequate bandwidth in order to accommodate the maximum 
possible Doppler shift of the CPM wave. The maximum possible Doppler 
shift depends on the vehicle speed, relative to the satellite. It is calculate 


as follows: 


^fmax — /o 


VrncLX 


(7.47) 
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where /o is the carrier frequency, u„iax is the maximum vehicle speed, and c is 
the speed of light. As stated earlier, the 7-power device produces an apparent 
Doppler shift (as seen by NBFl and NBF2) equal to 7A/. Therefore, NBFl 
and NBF2 must be wide enough to accoimt for the maximum apparent 
Doppler shift, In Fig- 7.13, the system parameters are listed. From 

those parameters and equation (7.47), we calculate ^fmax — 138 Hz, and 
7 A/max = 1104 Hz. Therefore, NBFl and NBF2 must allow for a frequency 
shift of ±1104 Hz, so they must have bandwidths of at least 2208 Hz. On 
Fig. 7.17 they have bandwidths of 3200 Hz, which is more than enough to 
allow for the Doppler shift, and narrow enough to reject the higher order 
sinusoids that appear at spacings of iZ* above the n 1 line. 

Since NBFl and NBF2 must be quite wide, their memory is reduced con- 
siderably, and they Me more susceptible to noise. This effect can be viewed 
as follows. Every time a symbol of {±1} is transmitted, the instantaneous 
CPM frequency (the derivative of the phase) is in the center of NBFl, so the 
filter gets excited. The filter then begins to resonate at its center frequency, 
with the filter output decaying exponentially after the data symbol ends. 
The filter continues to decay until it gets excited again, by another {±1} 
symbol transmission. Similarly, NBF2 gets excited when a symbol of {-1} 
is transmitted, and decays between {—1} s3rmbol tra n s m issions. The decay 
can be described by a time constant of the form 1/2x5, so a wide filter 
has a short time constant, and decays quicker than a narrow filter. Once 
the filters are completely decayed, the phase information they contain is no 
longer useful, so the recovered phase tends to drift randomly between 0 and 
2x. In order to limit these effects, we can narrow the post-detection filter, so 
that it has a long time constant, and retains the phase information longer. 

Viewing the action of the filters in the frequency dom^, we can visu- 
alize the CPM signal as a random process that contains spectral lines and 
continuum. Filters NBFl and NBF2 select the appropriate spectral lines, 
while allowing some amoimt of continuum to pass through. This continuum 
is then a source of noise, in the sense that it is an undesirable signal. If the 
filters are narrow, the signal that passes through them is predominately si- 
nusoidal, but if the filters Me wide, more noise passes through. The amount 
of spectral line power relative to continuum at the center frequency of filter 
NBFl is then related to the frequency of {±1} symbol transmissions, pulse 
shape, M, and h. SimilMly, the power of the spectral line at the center of 
filter NBF2 is related to the frequency of {-1} symbol transmissions. The 
final filter, from this viewpoint, limits the amount of noise that passes to 
the synchronization signal. 
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The other mo6xZe-related effect that is magnified by the presence of the 
7-power device is the dynamically-changing fading. First, the magnitude of 
the fading is raised to the 7 power, causing a large dynamic range expansion. 
Second, the fading-induced phase distortions are multiplied by 7, increasing 
the swing of the phase. Figure 4.3 is the spectrum of the CPM signal after 
the 7-power device for a 1-RC pulse shape, when noise (10 dB SNR.) and 
fading {K = 10 dB) are added to the CPM signal. Notice the asymmetry 
of the spectrum; most importantly, where we used to see only spectral lines, 
we now see that the lines are modiilated by the fading. Therefore, when we 
recover the synchronization signal, we get a sinusoid that is phase-modulated 
by the fading channel. Therefore, the recovered phase contains a lot of 
jitter. Figure 4.4 shows the recovered phase under the above mentioned 
conditions. It has an initial transient, which seems strange, coming out of 
such a narrow filter. However, the initial transient can be explained by the 
fact that the final filter is slow to respond to an input signal, so its output 
is slow to increase in magnitude. In the complex plane, that corresponds 
to an output which is near the origin. Being near the origin, the phase can 
easily move from one quadrant to another, causing large swings, until the 
output magnitude increases. Once the magnitude of the filter output builds, 
its output phase is more stable. Although the mean value of the phase is 
correct (0 radians), the variam.ce is quite large, so the tracking is not very 
good. Shortly, we will describe a device that linuts these effects. 

Finally, the non-linearity increases the effects of noise in the original 
CPM signal by producing additional cross-terms between the noise and the 
desired signal. This tends to degrade the performance of the system. Once 
again, a narrow post -detection filter will limit these effects, and wiU be 
discussed in the next section. 

7.3.3 Doppler Shift 

As mentioned in the previous section, Doppler shift is an important factor in 
a mobile environment, especially when it is time- varying, due to the ch^ge 
in the relative motion between the vehicle and the satellite. However, since 
Doppler shift affects aU of the spectral lines equaUy, the timing info^tion 
is not lost. In other words, when the two Doppler-shifted spectral lines are 
multiplied together and the symbol timing is recovered, the shift is self- 
canceling. 

Because filters NBFl and NBF2 must be wide, however, to allow the 
CPM signal to shift by Af^ax, they decay quickly, after being excited. 
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R, = 2400 Hz 

T, = 1/2400 sec 
h= 1/8 

7 = 8 

M = 8 

g{t) = 1-RC 

ai € {±1, ±3,±5,±7} 

/o = 1.5 GHz 

v^ax = 100 km/hr 


Figure 7.13: System Parameters 


Therefore, the recovered phase is prone to drifting, because the memory 
in the filters is relatively short, and once the filters are completely decayed, 
the phase they produce is no longer reliable. This problem was overcome by 
using a very narrow final post-detection filter, which liimts the excursions in 
the recovered phase that pass through the relatively wide NBFl and NBF2. 


7.3.4 Reducing fading effects 

Basically, the fading is caused by reflections and blockages of the CPM signal 
due to obstructions in the sign 2 d path (i.e. buildings). Also, because the sig- 
nal path keeps changing as the vehicle moves, the received composite sipial 
changes over time. This makes it difficult to maintain a lock on the received 
signal. A step that greatly reduces the fading effects is to insert a hard lim- 
iter at the front end of the synchronizer. This device maintains a constant 
signal amplitude, but it produces additional frequencies at multiples of the 
carrier, while keeping the original signal intact. The additional frequencies 
must then be filtered out through BPFl. The hard-limiter thus limits the 
effect of fading to only phase- shifting of the received signal, rather than 
phase and magnitude distortion. In effect, it restores the constant-envelope 
characteristic of the CPM wave, leaving residual phase modulation. 

7.3.5 System Parameters 

The parameters that were used in the simulation, described both in the 
2 inalog emd digital domain, are shown in Fig. 7.13. 
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Parameter 

Analog 

Digital 

Sampling Frequency 

/, = 38.4 kHz 

/, = 16 

Symbol Time 

T, = 1/2400 sec 

T, = 1 

Doppler Shift 

= 138 Hz 

Afmax = 0.0579 

NBFl Bamdwidth 

3200 Hz 

4/3 

NBF2 Bandwidth 

3200 Hz 

4/3 

LBPF Bandwidth 

50 Hz 

1/48 

BPFl Bandwidth 

5000 Hz 

50/24 


Figure 7.14: Simulation Parameters 

The modified block diagram of the synchronizer is shown in Fig. 7.17. 
Most notably, it cont 2 ims the awided hard limiter, the narrow post-detection 
filter, and eliminates the divide-by-n, because n is always 1. In the simula- 
tion, the addition of the hard limiter was handled very simply. The polar 
form of the received CPM signal is: 

y{t) = R{t) cos[2x/of + 0(t)]. (7.48) 

The limiter simply keeps the magnitude of this signal constant, setting R{t) 
equal to 1. Therefore, (7.48) simplifies to 

y[t) = cos[2x/of -t- ^(i)]' (7.49) 

The next section will present the results obtained in simulations of the com- 
plete modified system, as shown in Fig. 7.17. 
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7.4 Results and Conclusions 

In [4] an analysis was presented that calculated the timing jitter of the 
recovered symbol timing, using the block diagram in Fig. 7.1. The main 
assumption in the analysis was that the spectnim of the integer-h CPM 
(after the 7-power device) consisted of a spectral line in the presence of 
a continuum, md that the continuum was constant near the spectral line. 
Therefore, the analysis calculated the magnitude of the spectral line, and 
the Vcdue of the continuum at the frequency of the spectral line. Then 
it showed how to calculate the timing jitter from those two vjJues. The 
above assumptioTi is very good, provided that the bandwidths of NBFl and 
NBF2 are small, compared to the symbol rate. Some results ^lre tabulated 
in Fig. 7.18 using a bandwidth of 2.0/T, for BPFl, and a bandwidth of 
0.005/Ti for NBFl and NBF2. Note that the jitter variance is normalized 
to (2x)^rad^, which means that when we calculate the standard deviation 
of the phase jitter, by taking the square root of the variance, we should 
interpret the units ats being fractions of a symbol time. 

7.4.1 Performance Criteria 

With that interpretation of phase jitter in mind, we need to establish ein 
acceptable level of jitter, in order to judge the performance of the synchro- 
nizer. Ultimately, we eire concerned with a method of producing reliable 
communication using this digital system, so we need a low bit-error rate. 
The relationship between error rate and phase jitter takes the form; 

Pe(At) = j P{e\At)f{At)dAt (7.50) 

where At represents the error in the timing, and f{At) is the probability 
distribution of At. A higher phase-jitter variance yields a higher timing- 
jitter variatnce, which means that (7.50) will produce a higher bit-error rate. 
Generally, a phase jitter corresponding to about 5% of a symbol time is 
considered good, and 10% of a symbol time is acceptable. Therefore, our 
goal is to keep the timing jitter below 10%, which means that the normalized 
phase jitter should be below (0.1)^ = 0.01. 

For narrow filters (NBFl and NBF2), the final post-detection filter need 
only serve as a zonal filter, allowing frequencies in the vicinity of Rs to 
pass. Therefore its bandwidth is of little concern. Fig. 7.19 presents the 
results when those filters eire much wider. These results were calculated 
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BPFl Hard. Limiter BPF2 



pag. 58, Report jf: 5 


Figure 7.17: Modified Symbol Synchronizer 
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Pulse Shape 

Ek/No (dB) 

Normalized Variance 

1-REC 

0 

unreliable 


10 

1.0 


20 

4.5 X 10-^ 


30 

3.1 X 10-“* 


40 

2.9 X 10-=* 


50 

2.7 X lO-*^ 

1-RC 

0 

unreliable 


10 

8.8 X lO--* 


20 

3.3 X 10-“* 


30 

2.1 X lO-"* 


40 

2.0 X 10-^ 


50 

1.8 X 10“' 


Figure 7.18: Analytical Results (NBFl and NBF2 bandwidth = 0.005/Ti) 
{M = S,h= 1/8, n = 1) 


using bandwidths of 2.0/7'j for BPFl, and a bandwidth of Q.S/Tf, for NBFl 
and NBF2. 

Notice that, for very large bandwidths of NBFl and NBF2, the calcu- 
lated jitter variance is unreliable because it falls out of the range 0 to 1, 
corresponding to phase jitter in the range 0 to 2t. The reason for that 
behavior is that the original assumptions which supported that particular 
analysis fall apeirt (i.e., the assumption that the continuum is constant near 
the spectrad line). In this situation, the iinal post- detection filter serves a 
far more important role, emd must be narrowed considerably. 

7.4.2 Pulse shape selection 

Even with a narrow post-detection filter, the system did not perform reliably 
with some choices of pulse shape, even without noise or fading! An example 
is with the 1-REC p\ilse shape. A typical graph of the recovered phase in 
this situation is shown in Fig. 7.21 without noise or fading. In this situation 
the phase drifts slowly from — ir to ir, without tiny sign of stabilizing. Even 
though it seems like there are sudden phase discontinuities at several points 
in the graph, these discontinuities are just shifts between — ir and t, which 
really cire not shifts at eill; they occur due to the definition of the arctangent 
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Pulse Shape 

Eb/No (dB) 

Normalized Vciriance 

l-REC 

0 

unreliable 


10 

unreliable 


20 

4.5 X 10~‘ 


30 

3.1 X 10“^ 


40 

2.9 X 10-^ 


50 

2.7 X 10-“* 

1-RC 

0 

imreliable 


10 

unreliable 


20 

3.3 X 10“^ 


30 

2.1 X 10-^ 


40 

2.0 X 10-“ 


50 

1.8 X 10-^ 


Figure 7.19: Analyticed Results (NBFl and NBF2 bandwidth = 0.500/T^) 
(M = 8, h = 1/8, n = 1) 


function. When the drifting phase gets close to one of the extremes and 
crosses over into the next quadrant in the x-y plane, the arctangent function 
redefines the signal to fit the new quadrant. In any case, the behavior of 
this waveform is unacceptable. 

The reason for the unacceptable behavior is due to the spectrum of 
the integer-h CPM signal, given the 1-REC pulse. We will see, in the next 
section, an analysis that calculates the spectrum, but basically it shows that 
the power in the spectrzd line at n = 1 is not strong enough to track reliably, 
given that the wide filters pass a Icirge amount of continuum. If the filters 
were narrow, then there would be no problem with the REC waveform, but 
because they are wide, they pass a lot of power that is not at the frequency 
of the spectral line of interest. Therefore, we need to select a ptilse that has 
a high power in the n = 1 spectral line. 

Power Spectrum Calculated 

In [9] Aulin and Sundberg present a method for calculating the spectrum of a 
CPM signal. First, they derived an exact expression for the autocorrelation 
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ftmction of the CPM signal: 

iZ(r) = R{T' + mT,) (7.51) 

- -T Tn+l iVf — 1 , 

_ f ‘ ^ p^^2irhk[q{t+T’-{n-m)T,)-g(t-nT,)]y^ 

-^0 n=l-t k=-{M-l) 

(7.52) 

over the interval: 

0 < T < (I + 1)T, (7.53) 

where pk is the probability of the symbol occurring. Note that the time 
difference t has been written as: 

T = r' + mT,; 0<r'<l, m = 0,1,... (7.54) 

The power spectrum is then the Fourier transform of the autocorrelation 
fimction: 

5(/) = 2Re{/ iZ(r)c--'2"7rj^ 

p—j2irfLTs r 'o / 

M-1 

c. = 

k = -(Af- 1) 
k odd 

In this partictilar situation, h is equal to 1 (after the power-law device). 
This yields \Ca\ = 1. It is shown in [9] that this situation causes the auto- 
correlation function to be periodic outside |r| = TT,. Therefore, R{t') can 
be written as: 

R{t) = Rcon{r) + Rdis{r) (7.57) 

where Rdis{T) is periodic, causing the discrete components in 5(/). In [12] it 
is also shown that when \Ca\ = 1, Rcon{r) = 0; |r| > LT,. In this instance, 
Rco7x{t) yields the continuous part of the spectrum, and Rd,a(r) yields the 
discrete components. Therefore, the complete spectrum is found by taking 
the Fourier transform of Rcoti{t), and finding the coefficients of the Fourier 
series of f?i,j(r). 
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To determine 'ise the fact that ^cotiCt) = 0; [rj > LTj. 

This means that iZd„(r) = iJ(T); |t| > XT,. It also means that Rdia(r) 
is completely defined over a period of 2T„ since we already know that the 
spectrum has discrete components at Rsj2. Therefore, we calculate 
over the range 2LT, < r < 4XT„ and extend it in either direction to get the 
complete signal. Now, subtracting Rdui’’’) from R{r) gives Rconir )• 

The continuum part of the spectrum is calculated by taking the Fourier 
transform of Rconi'’’)’ 

Sconif) — f Rconi’’’)^ ^ ^ 

J — oo 

= 2Re|j^^^iic 

where the limits of integration are reduced due to the known characteristics 
of Rcon ( ^ ) • 

Finally, the Fourier series coefficients of calculated: 


on(r)e-^2"-^"dr| (7.58) 


cr. = ^ Rdis{r)e-^^^^-^dr (7.59) 

^ Jo 

The results of numericedly evaluating equations (7.58) and (7.59) are plotted 
in Figs. 7.22 and 7.23, with the magnitude of the spectral lines equal to the 
power in each one. 

The spectrum for the 1-REC pulse is shown in Fig. 7.22, and for the 
1-RC pulse in Fig. 7.23. Notice how much greater the ratio of the spectral 
line to continuum is in the 1-RC case than in the 1-REC case. Using 1- 
RC gives a 9 dB gain in this ratio. The ideal case (ideal in the sense that 
it produces the strongest spectral lines) is an impulse, which produces a 
differential phase shift keyed signal (DPSK). In that case, 


«(i) = 


When integrated, g(t) becomes: 


g(f) = 


0; 0 < i < ^ 
i. ^ < t < T 


(7.60) 


(7.61) 


In this system, the ph£ise jumps abruptly at T,/2, cuid always by ^ln odd 
integral multiple of t. So, 5.) eissumes values that always increment by 
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TT over a symbol time, producing a signal whose complex envelope is a square 
wave. Therefore, the DPSK spectrum consists only of spectral lines, with 
no continuum, since a square wave can be expressed as a sum of sinusoids. 
Given the above, an intuitive reason why 1-RC produces a better spectrum 
(as far as tracking performance goes) than 1-REC is that y(t) is closer to an 
impulse for the 1-RC pulse than for the 1-REC pulse. 

However, this change in pulse shape is not without problems. The orig- 
inal CPM spectrum using 1-RC has a wider main lobe than with 1-REC. 
The sidelobes, however, aie eventually lower, partially reducing the affect 
of the wider main lobe on the out-of-band power. This seems to be a minor 
price to pay, given the ease of synchronization of the recovered signal. 


Filter action 

Given the relatively strong spectral lines of the 1-RC system, we can now 
illustrate the action of NBFl and NBF2. As referenced in earlier chapters, 
NBFl gets excited periodically, then decays between excitations. For M = S, 
it gets excited an average of 1 time in 8, specifically when a {-1-1} symbol 
is transmitted. This phenomenon is illustrated in Fig. 7.24 for NBFl, and 
in Fig. 7.25 for NBF2, using one-sided filter bandwidths of 2R,/3 = 1600 
Hz for each filter. In the ideal situation, we expect each filter to oscillate 
continuously at its center frequency, which means that the complex envelope, 
relative to its center frequency, should remain constant. Obviously, from 
the graphs, this is not happening. In our situation, the filters decay greatly 
between excitations, and especicdly between symbols. Because of the short 
decay time, the filters lose their memory quickly, and the phase is no longer 
useful after the filters are completely decayed. At that point, the phase 
begins to drift, but adding a narrow post-detection filter alleviates some of 
the phase drifting by increasing the memory of the effective combination of 
filters NBFl, NBF2, amd LBPF. 

It is fairly obvious that a smaller alphabet size, Af, would increase the 
likelihood that NBFl would be excited (a {-1-1} symbol would be sent more 
frequently), eiUowing a shorter decay time between excitations. This is equiv- 
alent to increasing the spectral line power in the n = 1 line. This would 
improve the tracking performance also, but at a cost of lower information 
throughput. 
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Pulse Shape 

Eb/No (dB) 

K (dB) 

AF (Hz) 

Normalized Jitter Variance 

1-RC 


60 

0 

1.0 X 10--^ 


5 

60 

0 

1.8 X lO-** 


10 

60 

0 

6.9 X 10-^ 


120 

10 

0 

1.6 X lO-^* 


120 

13 

0 

1.3 X 10-*^ 


10 

10 

0 

1.1 X 10-“ 


10 

10 

150 

1.3 X 10"“ 


Figure 7.20: Final simulation results 


Recovered phase 

Figs. 7.26 and 7.27 show two examples of phase recovery, the only difference 
being in the delay time between when the signal is transmitted and when it 
is received. In Fig. 7.26, a zero delay is used, producing a recovered phase 
whose average value is 0. In Fig. 7.27, a delay of 0.25Tj is used, producing a 
phase whose average value is — t/ 2. Note that a delay of Tj is equivalent to 
zero delay, and would produce a recovered phetse of 0. Since the synchronizer 
produces an output that is proportional to the delay time, it is now clear that 
the synchronizer is performing its intended task. In order to illustrate the 
effects of Doppler shift, noise magnitude, and fading severity on the variance 
of the recovered phase, some s imula tion results are tabulated below. 

Interestingly, the variance for the no-noise, no-fading, no-Doppler shift 
case is not exactly zero, though it is quite small. The reason is due to data 
nois€y which stirfaces because of the randomness of the transmitted data 
symbols. Most importantly, it is related to the probability and of {+1} 
and {—1} symbol transmissions, and the autocorrelation function of the 
symbol generating source function (usually the source symbols are assumed 
to be xmcorrelated with one another.) A deterministic case, in which the 
symbol pattern is periodic (alternating {+!}, {^1} symbols, for example) 
does produce zero variance, because the filters (NBFl and NBF2) reach a 
steady-state after a few symbols are sent. 

In the worst-case noise situation without fading, Ef,/No = 0 dB, and 
the jitter variance is 10“^. This means that, on average, the timing jitters 
about ±(10“^)^/^ = ±0.032 of a symbol time, which is acceptable. In the 
worst case fading situation, without noise, if = 10 dB, and the timing jitters 
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about ±(1.6 X 10"®)^/^ = ±0.0013 of a symbol time. This suggests that, in 
the expected environment, the noise is more of a problem than the fading. 

Several combinations of noise and fading levels are included in Fig. 7.20, 
but perhaps the most Lmportaint two are the ones closest to the expected 
case; namely the situation when EfffNo = 10 dB, JiT = 10 dB, with varying 
values of A/. Without any Doppler shift (A/ = 0), the variance of the 
phaise jitter is 1.1 x 10““*, which means that the timing typically jitters 
±(1.1 X 10““*)*/^ = ±0.01 of a symbol time, which is acceptable. With a 
maximum Doppler shift of 150 Hz (larger than we expect), the timing jitter 
increases sUghtly, to ±(1.3 X 10-'*)*/^ = ±0.011 of a symbol time. The reason 
why the value is slightly higher for the Doppler-shifted case is a combination 
of several factors. First, the spectral lines are shifted close to the edge of the 
filters (NBFl and NBF2), so their amplitudes are starting to be reduced, due 
to the rolloff of the filters at the edge of their bandwidths. Also, the filters 
used in the simulation are HR, which means they do not have linear phase, 
but they are easier to use, in terms of memory requirements and speed of 
execution. Since they are not linear phase, especially at the outer edges 
of the frequency band, the signals that are near the edge experience some 
additional phase distortion, which, combined with the amplitude reduction 
at the filter edge, increases the jitter variance. In addition, a factor that is 
always present is simulation ’’noise” which translates into slight uncertainty 
in the accuracy of the simulation results, due mostly to finite averaging time 
in the simulation. However, even with the 11 percent increase in standard 
deviation between the no-Doppler and maximum-Doppler cases, the system 
is still quite useful. 
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1-RC 
M = 8 
h = 1/8 

Eb/Nfl = 10 dB 
K= lOdB 






Chapter 8 

Simulation results 


In this Chapter we present simulation results that allow a comparison be- 
tween different classes of trellis-encoded continuous- phase modulation (CPM) 
schemes with differential demodulation and optimum (maximum-likelihood, 
i.e., Viterbi) decoding. The performance analysis will be based on bit error 
probability versus the ratio SiyjNo between the average energy per informa- 
tion bit and the noise power spectral density of the additive white Gaussian 
noise aifecting the transmission. The baseline for our comparison has been 
chosen to be the performance of differentially-demodulated phase-shift key- 
ing (DPSK) with 8 levels. In fact, DPSK offers constant-envelope signals 
(like CPM) and is demodulated differentially (like our scheme). DPSK turns 
out to be more power-efficient than CPM, but on the other hand it occupies 
a broader bandwidth than CPM: thus, this comparison will show how our 
trellis- encoded CPM scheme will trade bandwidth efficiency for power effi- 
ciency. However, as will be demonstrated, trellis coded CPM comes fairly 
close, and indeed somewhat better, in performance to the filtered DPSK par- 
ticularly in a iZictan fading environment. 

Out simulation will include the effects of additive white Gaussian noise 
chaimel, Rician fading, and Doppler frequency shift. Our goal is to get the 
best combination of 

o Frequency shaping pulse g{t). 

o Full- or partial-response CPM. 

o Receiver filter shape and bandwidth. 

o Trellis- coded modulation (TCM) scheme. 
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Hereafter, the modxilation index h for CPM has always been chosen equal 
to 1/M, M the number of symbol carried by CPM. We hasten to observe 
that if the demodulation were optimum, i.e., coherent and based on Viterbi 
algorithm, then the optimum choice of h would involve a subtle tradeoff 
between power eind bandwidth. With differential detection, it is known that 
increasing h will increase the phase difference among CPM signals associated 
with different symbols, and hence increase the signal distance, which will 
ultimately restilt into a smaller error probability. Since increasing h will 
adso increase the bandwidth occupancy, the strategy that we have chosen 
is the following: choose h = 1/M, and the phase pulse shapes that satisfy 
the bandwidth constraints for that value of the modulation index. Soft 
decision demodulation was optimized based on analysis of the eye patterns. 
A complete set of precoded partial response signals have been studied under 
the AWGN model. It was observed that partial response systems shows a 
performance degradation larger than full-response CPM. For this reason, 
full-response CPM was analyzed more thoroughly with the aim of choosing 
the final candidate system. The selected full-response scheme was then 
tested \mder the Hician fading environment. 


8.1 Description of simulation results 

In this section we present the best combinations of TCM schemes, continuous- 
phase modulations, and receiving filters for a number of channels. The fol- 
lowing situations are considered: 

1. Frequency pulse shapes as shown in Table 8.1. In particular: 

• LREC denotes rectangular frequency pulse with duration LT, (in 
particular, L = 1 denotes full-response signaling). 

• LHCS denotes “half-cycle sinusoidal” frequency pulse with dura- 
tion LTj. 

• LSRC denotes “spectral raised cosine” with duration and 
“rolloff” factor a. 

• GMSK denotes a Gaussian frequency pulse. This has a theoreti- 
cad infinite duration, and is truncated to duration XT,. 

2. Receiver filter shapes: 

• LREC filter. 
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• Butterworth filter (unless otherwise specified, the Butterworth 
considered here has order 6 and a 3-dB bandwidth of 4800 Hz). 

• LaRC (raised cosine) filter (see Table 8.2). 

• Lb SRC (spectral raised cosine) filter (see Table 8.2). Here fij is 
the roUofF factor, zind oq is a normalization constant chosen such 
that //i2(f)df = 1. 

• Gaussian filter 

3. TreUis-coded modiilation schemes: 

• Rate-2/3, 4 states (see Fig. 1) 

• Rate-2/3, 8 states (see Fig. 2) 

• Rate-2/3, 16 states (see Fig. 3) 

4. Interleaving/deinterleaving sizes and depths: 

• No interleaving 

• (128,16) interleaving 

• (256,16) interleaving 

It should be observed here that we limited our consideration of TCM to 
schemes with no more than 16 states, and this for two reasons: 

1. While it is known that by increasing the number of states the perfor- 
mance of TCM increases, it has been observed that above a certain 
number of states the returns are diminishing, thus rendering less and 
less attr 2 ictive the increase in complexity entailed. 

2. We are interested in system performance at relatively small values of 
signal-to-noise ratio (typically, they correspond to error probabilities 
aroimd 10“^. In this range, TCM schemes with a large number of 
states perform more poorly than simpler schemes. In fact, the higher 
is the number of states, the higher is the signal-to-noise ratio necessary 
to take full advantage of the potential coding gain generated by the 
introduction of trellis- coded modulation. 

Interleaving the symbols eifter the TCM encoder and deinterleaving them 
before the TCM decoder would make the equivalent channel, as seen by 
the TCM scheme, more similar to a memoryless one. Since the decoding 
strategy is based on the simplifying assumption that the channel is indeed 
memoryless, it is expected that interleaving prove beneficial in the present 
context. 
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Table 8.1: Frequency pulses used in the simulation. 



Table 8.2: Impulse responses of filters pulses used in the simulation. 
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8.1.1 Gaussian channel 

In this section simulation results will be provided in the form of error prob- 
ability chauts corresponding to a transmission system operating over the 
additive white Gaussian noise (AWGN) chaimel, with intersymbol interfer- 
ence generated by non-ideal filters. 

• Fig. 8.4 shows the performance of 4-state, TCM-encoded CPM signals 
with full response signaling, HCSl frequency pulse, Af = 8 levels, h = 

0.125, and differential detection. This scheme carries 2 information 
bits per symbol. Here we compare the effect of different receiver filters, 
namely: 

1. IRC 

2. 2RC 

3. 2RC with (256,16) block interleaving 

4. SRC with Lb = l-Ol and rolloff = 0.4 

5. Same as above, with (256,16) block interleaving 

6. Butterworth 

7. Same &s above, with (256,16) interleaving 

Comparison of the performance with and without interleaving shows 
the benefit of introducing the latter. SRC filters provide the best 
performance in this situation. 

• Fig. 8.5 shows the performance of 4-state, TCM-encoded CPM signals 
with partied response signaling, 2REC frequency pulse, Af = 8 levels, 
h = 0.125, and differential detection. This scheme carries 2 informa- 
tion bits per symbol. Here we compare the effect of different receiver 
filters, namely: 

1. IREC 

2. IRC 

3. SRC 

4. Gaussian filter with 3-dB bandwidth 0.95 

5. SRC with Lb = 1-01 and roUoff 0j = OA 

This chcirt shows that the 2REC pulse offers poor performance for this 
system. 
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• Fie. 8.6 refers to the same situation as in Fig. 8.5, with the only 
change of an 8-state trellis in lieu of a 4-state trellis. The mcrease 
in complexity does not help. 

• Fie. 8.7 refers to the same situation as in Fig. 8.6, with the complexity 
of the TCM scheme further increased to 16 states. The same conclu- 
sion reached before about 2REC is still vaUd under the new conditions. 

• Fig. 8.8 shows the performance of 4-state, TCM-encoded CPM sign^ 
with partial response signaling, 2RC frequency pulse, M = 8 levels, 
h = 0.125, and differential detection. This scheme carries 2 informa- 
tion bits per symbol. Here we compare the effect of different receiver 

filters, namely; 

1. IREC 

2. IRC 

3. 3RC 

4. Gaussian filter with 3-dB bandwidth 0.95 

5. SRC with Xfl = 1.01 and rolloff /?/ = 0.4 

This chart shows that the 2RC pulse offers a performance for tins 
system which is better than 2REC, but still worse than the HCSl 
pulse considered before. 

• Fig 8.9 refers to the saime situation as in Fig. 8.8, with the only 
change of an 8-state treUis in Ueu of a 4-state trellis. The increase 
in complexity causes only a marginal performance improvement. 

• Fig. 8.10 refers to the same situation as in Fig. 8.8, with the complexity 
of the TCM scheme further increased to 16 states. The performance 
improvement is modest. 

• Fig. 8.11 shows the performance of 4-state, TCM-encoded CPM signals 
with partial response signaling, GMSK frequency pulse with BT = 0.4 
and truncated at LT, = 2, M = S levels, h = 0.125, and differential 
detection. This scheme carries 2 information bits per symbol. Here 
we compare the effect of different receiver filters, namely: 

1. IREC 

2. IRC 
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3. 3RC 

4. Gaussian filter with 3-dB bandwidth 0.95 

5. SRC with Lb = 1-01 and rolloff /3/ = 0.4 

This chart shows that the GMSK pulse offers a performance similar 
to 2RC. 

• Fig. 8.12 refers to the same situation as in Fig. 8.11, with the only 
change of an 8- state trellis in lieu of a 4-state trellis. The increase in 
complexity causes no performance improvement. 

• Fig. 8.13 refers to the same situation as in Fig. 8.11, with the complex- 
ity of the TCM scheme further increased to 16 states. The performance 
improvement is modest. 

• Fig. 8.14 shows the performance of 4-state, TCM-encoded CPM sig- 
nads with partial response signaling, HCS2 frequency pulse, M = 8 
levels, h = 0.125, and differential detection. This scheme carries 2 
information bits per symbol. Here we compare the effect of different 
receiver filters, namely: 

1. IREC 

2. IRC 

3. 3RC 

4. Gaussian filter with 3-dB bandwidth 0.95 

5. SRC with Xb = 1.01 and rolloff /3/ = 0.4 

This chart shows that the GMSK pulse offers a performance similar 
to 2RC. 

• Fig. 8.15 refers to the same situation as in Fig. 8.14, with the only 

change of 8- state trellis in lieu of a 4- state trelhs. The increase in 

complexity causes no performance improvement. 

• Fig. 8.16 refers to the same situation as in Fig. 8.15, with the com- 
plexity of the TCM scheme further increased to 16 states. There is no 
evident performance improvement. 

• Fig. 8.17 compares the performance of two TCM schemes with differ- 
ent number of states (and consequently, different complexity) under 
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the same conditions, namely, TCM-encoded CPM signals with partial 
response signaling, 2RC frequency pulse, M = 8 levels, h = 0.125, and 
dilFerential detection. The receiver filter is SRC with Lb = 1-01 and 
roUofF/3/ = 0.4. It is seen that, at least for large enough signal-to-noise 
ratios, increase in complexity involves a performance improvement. It 
should be noticed, however, that this improvement depends on the spe- 
cific receiver filter involved, and different filters give different returns, 
eis shown by comparison of Figs. 8 and 9. 

• Fig. 8.18 compares the performance of two TCM schemes with differ- 
ent number of states (and consequently, different complexity) under 
the same conditions, namely, TCM-encoded CPM signals with par- 
tial response signaling, GMSK frequency pulse truncated at XT, = 2 
and with BT = 0.4, M = 8 levels, h = 0.125, and differential detec- 
tion. The receiver filter is IRC. The same conclusions reached in the 
examination of Fig. 8.17 stUl hold true. 

• Fig. 8.19 compares the performance of three TCM schemes with dif- 
ferent number of states (and consequently, different complexity) under 
the same conditions, namely, TCM-encoded CPM signals with partial 
response signaling, HCS2 frequency pulse, M = 8 levels, h = 0.125, 
md differential detection. The receiver filter is Gaussian with BT, = 

0.95. As expected, a 16-state scheme provides a larger coding gain, but 
this is true only at higher values of signal-to-noise ratio. The overall 
performance is slightly worse than with 2RC pulses. 

• Fig. 8.20 shows the performance of 4-state, TCM-encoded CPM signals 
with full response signaling, IREC frequency pulse, Af = 8 levels, h = 
0.125, and differential detection. This scheme carries 2 information 
bits per symbol. Here we compare the effect of different receiver filters, 
namely: 

1. IREC 

2. IRC 

3. 2RC 

4. Gaussian filter with BT, = 0.95 

5. SRC with Lb — l-Ol and rolloff /3/ = 0.4 

6. Butterworth 


156 


CHAPTER 8. SIMULATION RESULTS 


Butterworth filters provide the best performance in this situation, but 
the system performance is worse than with full-response HCSl pulse 
(See Fig. 8.4 above). 

• Fig. 8.21 refers to the same situation as Fig. 8.20, with the only change 
of an 8-state TCM scheme instead of a 4-state scheme. 

• Fig. 8.22 refers to the same situation as Fig. 8.20, with the only change 
of an 16-state TCM scheme instead of a 4-state scheme. 

• Fig. 8.23 shows the performance of 4-state, TCM-encoded CPM signals 
with full response signaling, IRC frequency pulse, M = 8 levels, h = 

0.125, and differenti2d detection. This scheme carries 2 information 
bits per symbol. Here we compare the effect of different receiver filters, 
namely: 

1. IREC 

2. IRC 

3. 2RC 

4. Gaussian filter with BT, = 0.95 

5. SRC with Lb = 1*01 rolloff = 0.4 

6. Same as above, with (256,16) block interleaving 

7. Butterworth 

The SRC filter with interleaving provide the best performance in this 
situation. The resulting performance is only slightly worse than that 
of 8DPSK. 

• Fig. 8.24 refers to the same situation as the previous figure, but it 
emphasizes the effect of interleaving. It is seen that for an error prob- 
ability of 10“^ interleaving improves the signal-to-noise ratio of about 
0.8 dB. 

• Fig. 8.25 refers to the same situation as Fig. 8.23, with the only change 
of an 8-state TCM scheme instead of a 4-state scheme. Also, inter- 
leaved SRC filtering is not considered here. 

• Fig. 8.26 refers to the same situation as the previous figure, but it 
emphasizes the effect of interleaving. It is seen that for an error prob- 
ability of 10“^ interleaving improves the signal-to-noise ratio of about 
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0.8 dB, although RC filters do not provide as good a performance as . 
SRC filters. 

. Fig. 8.27 refers to the same situation as Fig. 8.25. 

effect of interleaving. It is seen that for an error ^ 

interleaving improves the signal-to-noise ratio of about 0.8 ^B. Here 
SRC filters are considered, wWch improves the performance with re- 

spect to RC filters. 

. Fig. 8.28 refers to the same situation as Fig. 8.25. It eniphasizes the 
effect of interleaving and of the interleaving size and depth. 

• Fig. 8.29 refers to the same situation as Fig. 8.27, with the receiving 
filter changed from SRC into Butterworth. This proves once agam e 
benefits of interleaving, but the performance is not as good as wi 

SRC filter. 

. Fig. 8.30 compares 4-state and 8-state TCM schemes in the same sit- 
uation as in Fig. 8.28. Once again, it is seen that ^ mcrease m the 
number of states improves the performance at higher values 
signal-to-noise ratio. 

. Kg. 8.31 shows tho performance of 16-state, TCM-encoded CPM sig- 
with foil response signaling, IRC frequency pulse, M = 8 
k = 0.125, and differential detection. This scheme c^es 2 infor^ 
tion hits per symbol. Here we compare the effect of different receiver 

filters, namely: 

1. IRC 

2. 2RC 

3. 2RC with (256,16) block interleaving 

4. SRC with Lb = 1-01 rolloff /3/ = 0.4 

5. Same as above, with (256,16) block interleaving 

6. Butterworth 

7. Butterworth, with (256,16) block interleaving 

The SRC filter with interleaving provide the best performance 
situation. The resulting performance is only slightly worse than a 

of 8DPSK. 
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• In the same situation as in the previous figure, Fig. 8.32 compares 
interleaved and non-interleaved transnussion with SRC filtering. Here 
interleaving provides a gain of 0.6 dB in signal-to-noise ratio. 

• In the same situation as in Fig. 8.31, Fig. 8.33 compares interleaved 
and non-interleaved transmission with Butterworth filtering. Once 
more, interleaving provides a gain of 0.6 dB in signal-to-noise ratio. 

• Fig. 8.34 shows the performance of 4-state, TCM-encoded CPM signals 
with fun response signaling, GMSKl frequency pulse with BT, = 0.7, 
M = 8 levels, h = 0.125, and differential detection. This scheme 
carries 2 information bits per symbol. Here we compare the effect of 
different receiver filters, namely: 

1. IREC 

2. IRC 

3. 2RC 

4. Gaussian filter with BT, = 0.95 

5. SRC with Lb = 1-01 and roUoff /3/ = 0.4 

6. Butterworth 

The Butterworth filter provides the best performance in this situation, 
although the resulting performance is several dB worse than that of 
8DPSK. 

• In the same situation as in the previous figure. Fig. 8.35 considers an 
8-state TCM scheme. 

• In the same situation as in Fig. 8.34, Fig. 8.36 considers a 16-state 
TCM scheme. 

• Fig. 8.37 compares 4-, 8-, and 16-state TCM scheme for full response 
signaling, GMSKl frequency pulse with BT, = 0.7, M = 8 levels, 
h = 0.125, and differential detection. The receiver filter here is But- 
terworth. The same conclusions listed before about the effect of an 
increzise in the number of states of the TCM scheme hold here: such 
increase is beneficial only at high-enoUgh signal-to-noise ratio. 

• Fig. 8.38 shows the performance of 4-state, TCM-encoded CPM signals 
with full response signaling, HCSl frequency pulse, M = 8 levels, h = 
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0.125, and differential detection. This scheme carries 2 information 
bits per symbol. Here we compare the effect of different receiver filters, 
namely: 

1. SRC with Xb = 1.01 and rolloff = 0.4 

2. Same as above, with (256,16) block interleaving 

3. Butterworth 

4. Butterworth with (256,16) block interleaving 

The SRC filter with interleaving provides the best performance in this 
situation. 

• In the same situation as in the previous figure. Fig. 8.39 compares 
interleaved and non-interleaved transmission with SRC filtering. Here 
interleaving provides a gain of 0.9 dB in signal- to-noise ratio. 

• In the same situation as in Fig. 8.38, Fig. 8.40 compares interleaved 
and non-interleaved transmission with Butterworth filtering. Here in- 
terleaving provides a gain of 0.8 dB in signal-to-noise ratio. 

• Fig. 8.41 shows the performance of 8-state, TCM-encoded CPM signals 
with full response signaling, HCSl frequency pulse, M = 8 levels, h = 

0.125, emd differential detection. This scheme carries 2 information 
bits per symbol. Here we compare the effect of different receiver filters, 
namely: 

1. IRC 

2. 2RC 

3. 2RC with (256,16) block interleaving 

4. Gaussiem filter with BT, = 0.95 

5. SRC with Lb = 1.01 and rolloff /3/ = 0.4 

6. Same as above, with (256,16) block interleaving 

7. Butterworth 

The SRC filter with interleaving provides the best perform^mce in this 
situation. 

• In the S 2 ime situation as in the previous figure. Fig. 8.42 compares 
interleaved and non-interleaved transmission with SRC filtering. Here 
interleaving provides a gain of 0.9 dB in sign«il- to-noise ratio. 
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• Fig. 8.43 shows the performance of 16-state, TCM-encoded CPM sig- 
nals with full response signaling, HCSl frequency pulse, Af = 8 levels, 
h = 0.125, and differential detection. This scheme carries 2 informa- 
tion bits per symbol. Here we compare the effect of different receiver 
filters, namely: 

1. IRC 

2. 2RC 

3. 2RC with (256,16) block interleaving 

4. SRC with Lb = 1-01 and rolloff /3f = 0.4 

5. Same as above, with (256,16) block interleaving 

6. Butterworth 

7. Butterworth with (256,16) block interleaving 

The SRC filter with interelaving provides the best performance in this 
situation. 

• In the same situation as in the previous figure. Fig. 8.44 compares 
interleaved and non-mterleaved transmission with SRC filtering. Here 
(256,16) block interleaving provides a gain of about 1 dB in signal-to- 
noise ratio, while (128,16) interleaving offers only 0.4 dB. 

• Fig. 8.45 refers to the same situation as the previous figure, with But- 
terworth receiving filter. Comparison of interleaving schemes of dif- 
ferent size shows that (256,16) block interleaving provides a gain of 
about 0.8 dB in signal-to-noise ratio, while (128,16) interleaving offers 
only 0.4 dB. 

• Fig. 8.46 compares the performance of 4-state and 8-state TCM-encoded 
CPM signals with full response signaling, HCSl frequency pube, M = 

8 levels, h = 0.125, and differential detection. This scheme carries 2 
information bits per symbol. The receiver filter is Butterworth. 

In conclusion, for the AWGN chsnnel model (i-e., in the Absence of {&d- 
ing) an 8-state TCM-encoded CPM signals with full response signaling, IRC 
frequency pulse, M = 8 levels, h = 0.125, and differential detection, in con- 
junction with an SRC receiver filter with Lb = 1-01 and rolloff 0/ = 0.4, and 
(256,16) block interleaving, offers the best compromise between complex- 
ity and performance in the signal-to-noise range and among the candidate 
schemes considered in this Report. 
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8.1.2 Fading channel 

We now consider a channel ^l^Fected by additive Gaussian noise, intersymbol 
interference due to filtering effects, Doppler frequency spread, and Rician 
fading. The following par^uIleters were selected: 

1. Rician parameter K: chosen to be 10 . 

2. Doppler frequency spread : chosen to be 20 and 40 Hz. 

In our simulations we have examined the behavior of the best schemes 
found for the AWGN channel, i.e., IRC and IHCS. The latter has a higher 
degradation (as shown in Fig. 8.44), but is spectrally more efficient than 

IRC. 

The results obtained are illustrated in the figures that follow. 

• Fig. 8.47 illustrates the performance of the full-response 4-state trel- 
lis coded CPM with HCS pulse shape, SRC filtering, two different 
interleaving sizes, and a Doppler spread of 20 Hz. 

• Fig. 8.48 shows the performance of the previous system with an 8-state 
code. 

• Fig. 8.49 refers to the same situation as the previous figure. The only 
change is in the number of TCM states, that is 16 here. Here, a 
comparison is made between the proposed CPM system (which has 
constant envelope) and the filtered 16-state 8DPSK (whose envelope 
is not consteint). 

With (256, 16) interleaving the CPM system offers the same perfor- 
mance as the non-constcint envelope filtered 8DPSK with (128, 16) 
interleaving. This will make the CPM relatively superior the 8DPSK 
system. 

With (128, 16) interleaving, our CPM system shows a degradation of 
only .5 dB with respect to 8DPSK. In the presence of significant chan- 
nel nonlinearities due to power amplifiers driven at or near saturation 
for better power efficiency, CPM will perform better than 8DPSK be- 
cause in the latter scheme envelope fluctuations will cause phase fluc- 
tuations, which in turn will be reflected into further performance degra- 
dation. 
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• Fig. 8.50 shows the performance of the full response 4-state trellis 
coded CPM signal, over the Rician fading channel with a Doppler 
spread of 40 Hz, with HCS pulse shape, SRC filter, and two different 
interleaving sizes. 

• Fig. 8.51 demonstrates the expected 2 dB loss, for a 4-state trellis 
coded CPM signal, due to the presence of fading, with and without 
interleaving. 

• Fig. 8.52 is the same as Fig. 8.50 with an 8-state TCM scheme. 

• Fig. 8.53 compares the performance of the 8-state trellis coded full re- 
sponse CPM signal over AWGN channel with Rician fading and 40 Hz 
Doppler spread, with and without interleaving. This system shows a 
performance degradation, due to fading, of only .9 dB. 

• Fig. 8.54 is the same as Fig. 8.50 with a 16-state code. 

• Fig. 8.55 compares the performance of full-response CPM signal with 
HCS pulse shape over AWGN channel with Rician fading and a 40 Hz 
Doppler spread for 4, 8, and 16-state codes. A performance improve- 
ment of about 2 dB caused by the increase in the number of code 
states from 4 to 16 is observed here. 

• Fig. 8.56 compares the performance, over two different Rician fading 
rharmpls with the Same interlaving scheme, of the following two sys- 
tems: 

8 DPSK : 16-statetrelliscoded, filtered, over Rician fading with 20 Hz 
Doppler spread, and (128, 16) interleaving. 

Diflferential CPM ; 16-state trellis coded, over Rician fading with 
40 Hz Doppler spread, and (128, 16) interleaving. 

• Fig. 8.57 illustrates the performances of 4-state trellis coded CPM 
signal, over the Ricicin fading chaimel, with IRC pulse shape for two 
different interleaving sizes. 

• Fig. 8.58 illustrates the performancs of 8-state trellis coded CPM sig- 
nal, over the Riciain fading chaimel, with IRC pulse shape for two 
different interleaving sizes. Comparison is drawn with the 16-state 
non-constemt envelope interleaved, with the same Doppler spread of 
20 Hz, filtered 8DPSK, indicating a loss of only 1 dB. 
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• Fig. 8.59 illustrates the performancs of 16-state trellis coded CPM 
signal, over the Rician fading channel, with IRC pulse shape for two 
different interleaving sizes. Comparison is drawn with the 16-state 
non-constcint envelope interleaved, with the same Doppler of 20 Hz., 
filtered 8DPSK, indicating a loss of only 1.1 dB. 

• Fig. 8.60 is the same as fig. 8.57 except, the Doppler spread is 40 Hz. 

• Fig.8 61 is the same as fig. 8.58 except, the Doppler spread is 40 Hz. 

• Fig. 8.62 compares the performance, over two different Rician fad- 
ing channels with the same interleaving scheme, of the following two 
systems: 

8DPSK : 16-state trellis coded, filtered, over Rician fading with 20 Hz 
Doppler spread, and (128, 16) interleaving. 

Differential CPM ; 16-state treUis coded, over Rician fading with 
40 Hz Doppler spread, and (128, 16) interleaving. 

• Fig. 8.63 compares the performance, over the same Rician fading chan- 
nel with 20 Hz. Doppler spread, of: 

8DPSK : 16-state trellis coded, filtered, and (128, 16) interleaving. 
Differential CPM j 8 and 16-state trellis coded, and (128, 16) in- 
terleaving. 

We observe that, under the conditions of our simulation, the CPM scheme 
offers the same performance as a filtered 8DPSK signal, but with the addi- 
tional positive feature of a constant envelope. 71118, in essence, will make 
the CPM a more attractive system. We also observe that, in the presence of 
a finite-size block interleaving, larger Doppler spreads may cause a smaller 
performance degradation than smaller spreads. This is apparently due to 
the fact that Doppler frequency spread causes a reduction of the length of the 
bursts, which aie not completely interspersed by a finite-depth interleaver. 
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Performance of 4-state trellis coded 8-level CPM signal with HCSl pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.4 
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Performance of 4-staie trellis coded 8-Ievel CPM signal with HCSl pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.5 
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Performance of 8-state trellis coded 8-level CPM signal with 2REC pulse 
shape, h=0. 125, and differential detection for various receiver filtering. 


Figure 8.6 
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Performance of 8-state trellis coded 8-level CPM signal with 2REC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.7 
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Performance of 4-state trellis coded 8-level CPM signal with 2RC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 


Figtire 8.8 
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Performance of 8-state trellis coded 8-level CPM signal with 2RC pulse 
shape, h=0.125, and differenual detection for various receiver filtering. 


Figure 8.9 
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Performance of 16-state trellis coded 8-level CPM signal with 2RC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.10 
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Performance of 4-state trellis coded 8-level CPM signal with GMSK2, BT-.4 
pulse shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.11 
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Performance of 8-state trellis coded 8-level CPM signal with GMSK2, BT=.4 
pulse shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.12 
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Performance of 16-state trellis coded 8-level CPM signal with GMSK2, BT=.4 
pulse shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.13 
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Performance of 4-state trellis coded 8-level CPM signal with HCS2 pulse 
shape, h=0.125, and differenual detection for various receiver filtering. 


Figure 8.14 
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Performance of 8-state trellis coded 8-lcvel CPM signal with HCS2 pulse 
shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.15 
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Performance of 16-state trellis coded 8-level CPM signal with HCS2 pulse 
shape, h=0.125, and differential detection for various receiver filtering. 


Figure 8.16 
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Performance of trellis coded 8-Ievel CPM signal with 2RC pulse 
shape, h=0.125, and differential detection for 4 and 8-state codes. 
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Performance of trellis coded 8-level CPM signal with GMSK2, BT=.4 pulse 
shape, h=0.125, and differential detection for 8 and 16-state codes. 


Figure 8.18 
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Performance of trellis coded 8-level CPM signal with HCS2 pulse 
shape, h=0.125, and differential detection for 4, 8, and 16-state codes. 


Figure 8.19 
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Performance of 4-state trellis coded 8-level CPM signal with IREC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.20 
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Performance of 8-state trellis coded 8-level CPM signal with IREC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 
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Performance of 16-state trellis coded 8-Ievel CPM signal with IREC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.22 
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Figure 8.23 


Probability of Bit Error 


CHAPTER 8. SIMULATION RESULTS 


187 



Performance of 4-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figtire 8.24 
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Eb/No(dB) 

Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 
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Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.26 



Probability of Bit Error 


CHAPTER S. SIMULATION RESULTS 190 



5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 
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Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.27 
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Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.28 
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Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 
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Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection for 4 and 8-state codes.Trellis 
coded 8-level DPSK is also shown for comparison. 


Figure 8.30 
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Performance of 16-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0,125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 
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Performance of 16-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 
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Performance of 16-state trelUs coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.33 
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Performance of 4-state treUis coded 8-IeveI CPM signal with GMSKl. BT=,7 
pulse shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-Ievel DPSK is also shown for comparison. 


Figure 8.34 
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Performance of 8-state trellis coded 8-leveI CPM signal with GMSKl, BT=.7 
pulse shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.35 
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Performance of 16-state trellis coded 8-level CPM signal with GMSKl, BT=.7 
pulse shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 
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Trems coded 8-level DPSK is also shown for companson. 


Figure 8.37 
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Performance of 4-state trellis coded 8-level CPM signal with HCS 1 pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.38 
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Performance of 4-state trellis coded 8-level CPM signal with HCSl pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.39 
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Performance of 4-state mdlis coded 8-le.el 

Shane h=0 125, and differeniial detecuon with and wi*out interlea g. 
TrelUs coded 8-level DPSK is also shown for comparison. 
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Performance of 8-state trellis coded 8-IeveI CPM signal with HCSl pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.41 
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T, r 8-state trellis coded 8-level CPM signal with HCSl pulse 

125, and differential detection with and widtout interleaving. 
Trems coded 8-level DPSK is also shown for companson. 


Figure 8.42 
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Eb/No(dB) 

Performance of 16-state trellis coded 8-leveI CPM signal with HCSl pulse 
shape, h=0.125, and differential detection for various receiver filtering. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.43 
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Performance of 8-state trelUs coded 8-level CPM sign^ with HCSl pulse 
shape, h=0.l25, and differential detection with and without interleavmg. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.44 
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Performance of 16-state trellis coded 8-level CPM signal with HCS 1 pulse 
shape, h=0.125, and differential detection with and without interleaving. 
Trellis coded 8-level DPSK is also shown for comparison. 


Figure 8.45 
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Performance of trellis coded 8-level CPM signal with HCSl pulse 
shape, h=0. 125, and differential detection for 4 and 8-state codes. 
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Performance of 4-state trellis coded 8-Ievel CPM signal with HCSl pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
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Figure 8.48 






Probability of Bit Error 


CHAPTER 8. SIMULATION RESULTS 


212 


I j 

y 



5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

E|/No(dB) 


Performance of 16-staie trellis coded 8-level CPM signal with HCSl pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 16-state trellis coded 8-DPSK is also shown for comparison. 


Figure 8.49 
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Performance of 4-state trellis coded 8-level CPM signal with HCS I pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 


Figure 8.50 
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Figure 8.51 
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Performance of 8-state trellis coded 8-Ievel CPM signal with HCS 1 pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 
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Performance of 8-state trellis coded 8-level CPM signal with HCSl pulse 
shape, h=0.125, and differential detection over different channels. 


Figure 8.53 
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Performance of 16-state trellis coded 8-level CPM signal with HCSl pulse 
shape. h=0.125, and differential detection over Rician fading and AGWN 
channel. 


Figure 8.54 
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Performance of trellis coded 8-leveI CPM signal with HCSl pulse 
shape, h = .125, and differential detection, over Rician fading and 
AWGN channel, for 4, 6, and 16-state codes. 


Figure 8.55 
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Comparison between the performance of 16-state trellis coded 8-Ievel CPM 
signal with HCSl pulse shape, h=0.125, with differential detection and 
16-state 8DPK, over the Rician fading and AGWN channel, for two 
different Duppler spreads. 


Figure 8.56 
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Performance of 4-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 


Figure 8.57 




Probability of Bit Error 


CHAPTER 8. SIMULATION RESULTS 


221 



5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 

E(/No(dB) 


Performance of S-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 16-statc trellis coded 8-DPSK is also shown for comparison. 


Figure 8.58 
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Performance of 16-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 16-state trellis coded 8-DPSK is also shown for comparison. 


Figure 8.59 
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Performance of 4-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection over the Rician fading and AGWN 
channel, 16-state trellis coded 8-DPSK is also shown for comparison. 


Figure 8.60 
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Performance of 8-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection over Rician fading and AGWN 
channel. 16-state trellis coded 8-DPSK is also shown for comparison. 
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Performance of 16-state trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection over the Rician fading and AGWN 
channel. 16-state trellis coded 8-DPSK is also shown for comparison. 


Figure 8.62 
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Performance of trellis coded 8-level CPM signal with IRC pulse 
shape, h=0.125, and differential detection , over Rician fading and AGWN 
channel, for 4 and 8-state codes. 16-state trellis coded 8-DPSK is also 
shown for comparison. 


Figure 8.63 






Chapter 9 

Conclusions 


We have considered a satellite-based mobile communication scheme based on 
continuous-phase modulated (CPM) signals used in conjunction with trellis- 
coded modulation (TCM). We have proved that the synergy between TCM, 
which improves error probability, and CPM signals, which provide constant 
envelope and low spectral occupancy, provides a satisfactory solution to the 
problem of transmitting on mobile satellite chaimels. 

CPM signals are differentially demodulated symbol-by-symboL By doing 
this, the power efficiency of CPM is not exploited, but its spectral properties 
are retained, and interleaving/deinterleaving is made possible. 

This transmission scheme was extensively simulated under several condi- 
tions. Several frequency pulse shapes, receiver filter shapes, TCM schemes, 
interleaving/deinterleaving sizes and depths, and fading channels were con- 
sidered. It has been foimd that, over the AWGN channel model (i.e., in the 
absence of fading) an 8-state TCM-encoded CPM signals with full response 
signaling, IRC frequency pulse, M = 8 levels, h — 0.125, and differential 
detection, in conjunction with an SRC receiver filter with Lb — 1.01 and 
rolloff ^/ = 0.4, and (256,16) block interleaving, offers the best compro- 
mise between complexity and performance in the signal- to- noise range and 
among the candidate schemes considered in this Report. For fading chan- 
nels, full-response CPM with IHCS pulses offers the same performance as a 
filtered 8DPSK signal, but with the additional positive feature of a constant 
envelope. This, in essence, will make the CPM a system more attractive 
than DPSK. 

In conclusion, our results show that, for a careful selection of system pa- 
rameters, on satellite-based land mobile communication channels a differentially- 
detected CPM scheme can offer an error performance which is essentially the 
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same as differentially coherent PSK. Since PSK does not use bandwidth in a 
very efficient way, it has to be band-pass filtered to meet the requirements of 
closely-spaced mobile-radio channelization. As a result, its envelope is not 
constant, and its performance would be degraded by power amplifiers oper- 
ated at or near saturation for better power efficiency. Since CPM with the 
parameters chosen is bandwidth efficient it does not require narrow filtering, 
and consequently offers constant or near-constant envelope, thus preventing 
a possible cause of serious performance degradation. 
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Appendix A 

The simulation package 


A.l Structure of the simulation package 

The computer simulation package developed for the analysis ^d design of 
the trellis-encoded CPM scheme includes the followmg subpackages. 

1. Input/Output blocks. They input the P^T**"* 

system performance and store them mto an assigned file. They ako 
store the output quantities and store them into files for graphical pre- 

sent at ion. 

2. TCM encoder and Viterbi decoder. 

3 CPM modulator and demodulator. These blocks generate ^erent 
pulse shapes according to the input parameters, and demodulate the 
received signal by using a differential or non-coherent demodulator. 

4. Channel simulators. Addition of Gaussian noise and fading with or 
without shadowing. 

5. Filter simulators. Modules are available to simulate the 
families: IREC, LRC, Gaussian, Chebychev, Butterworth, and LS 
Gaussian and LSRC use the frequency-domain exact form, the other 
filter models are based on approximations. 

6. Simulator of the Doppler-frequency shift estimation and removal. 

7. Simulator of the timing recovery circuit. 
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Figiire A.l: Flow chart of the simulation program. 
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Fig. A.l shows the flow chart of the package as used to smulate the 
behavior of the TCM-encoded CPM system. The number of CPM sym o 
levels generated by the transmitter may be 2, 4, ar 8. The mary 
quaternary case have been validated by compa^g our simulation resets 
with results previously pubUshed. The pulse shapes Aosen were LREC, 
LRC, and GMSK. The binary cases considered are without preco g, an 
the samples at the output of the differential detector were taken at m^ger 

multiples of signal symbol time, i.e., = nT,. For for ^ 

preco^g was used on the source data. The precodmg rules axe described 

in the following: 


Transmitter: 


Receiver: 


output of data generator 
output of precoder 

input of CPM modulator 

detector output 


The decoding rules are: 


Vn = -6 or 2 
y„ = -4 or 4 
yn = -2 or 6 
Vn = 0 


bn € { 0 , 1 , 2 , 3 } 

On=bn- On-l 

Un ^ {0, 1»2,3} 

Cn = 20n — 3 
Cn€{±l,±3} 

y„ = C„ + Cn-1 = 2(o„ + On-l) ' 6 
= 2[(5n - o„-i mod 4) + On-i] - 6 
e{0,±2,±4,±6} 


bn = 0 
bn = l 

bn = 2 

bn = 2 


The relations among bn,On-i, and y„ are summarized in Table A.l. 

A.1.1 Multipath fading simulator 

In our description of the channel used in this Report, we have Hsted three 
fading models, namely, the Rayleigh, Rice, and Loo u^dels Here we de- 
scribe the structure of the software written to simulate these channel models, 

and provide their validation based on simulation results. 

The commonly accepted signal propagation model for satelhte-based mo- 
bile communication is based upon consideration of three signal paths from 
the satellite to the mobile, namely, a line-of-sight or direct path, a specu ar 
path, and a scatter path or multipath. The specular ^d direct path axe 
combined to form the coherent component, while the diffuse path forms the 
noncoherent component. 
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bn 

an-1 

Vn 

0 

0 

-6 

0 

1 

2 

0 

2 

2 

0 

3 

2 

1 

0 

-4 

1 

1 

-4 

1 

2 

4 

1 

3 

4 

2 

0 

-2 

2 

1 

-2 

2 

2 

-2 

2 

3 

6 

3 

0 

0 

3 

1 

0 

3 

2 

0 

3 

3 

0 


Table A.l: Precoding table for quaternary CPM. 
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The latter component is Gaussian, while the coherent component, in 
terms of its complex envelope, is either a constant (Rayleigh model. Rice 
model) or is subjected to a log-normal transformation (Loo model). The 
log-normal transformation represents the effects of foliage attenuation, and 
is usually referred to as shadowing. 

The sum of a constant-envelope signal and a Gaussian process has an 
envelope modeled by a Ricim probability density function. In the absence 
of the constmt-envelope component, the probability density function is of 
the Rayleigh type. The functional block diagram for the Rayleigh fading 
simulator is shown in Fig. A.2. The Gaussian noise source generates a 
sequence of white Gaussian noise samples. These samples pass trough a 
shaping filter, and modulate the in-phase and quadrature parts of a radio- 
frequency carrier. The theoretical model for the shaping filter is depicted in 
Fig. A.3, md its frequency response can be expressed in the form 

mfl.l A [1 

”'"lo |/|>A 

where = v/A is the Doppler frequency spread corresponding to a vehicle 
speed V when the carrier wavelength is A, aind ^4 is a gain constant. The 
probability density function of the fading sample is 

fa(a) = 2ae”®^, a > 0. 

The phase of the Rayleigh sample has a uniform distribution between -tt 
and T. 

The functional block diagram for the Rician fading simulator is shown 
in Fig. A.4. The upper arm is essentially the same as for the Rayleigh 
model simulator- The lower arm generates a Doppler- shifted signal compo- 
nent, which represents the direct path propagation. The probability density 
function for the Rici^ul envelope is 

^(a) = 2a(l + (^2a^ii:(l -h IT)) , a > 0, 

where K is the ratio between the power in the direct path and the power in 
the multipath signal. The phase probability density function of the Rician 
fading samples is given by 

/@(^) = -I- ;^v/^cos0e-'^”""^[2 - erfc(v^cos<?)]. 

27T 2y 7T 
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Figure A.2: Simulator for the Rayleigh model. 
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Figure A. 4: Simulator for the Rice model. 
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Figtue A.5: Simulator for the Loo model. 
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an samples, K = 10. 
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The block diagiem of the tor the lower 

Fig. A.5. The tunt^tor *kdowing component. 

/ .. \2 ^ 


i 


< 00 . 


„ I, weU known that a random process p(.) whose prohahOlt, denslt, fonc- 
« “ "tag normal can he represented in the form 

y(t) = exp{-X’(i)}> 

/t 2 TKe bandwidth of 

rsLptallt" — 

that of the Rayleigh component. 

Validation faded samples is used in the 

The baseband signal samples can be ob- 

computer simulations. Th , signal. The shaping filter for the 

tained by forming the bj . tod-order digital Butterworth 

RAjlelgh ma|4tnde of the fceqnency response is de- 

filter ’rt* Jianiw^ f„ to frequency response of the 

Srdlrder'digitai Butterworth SL^d^eTnot tre“ase 

noise components to pass through the Me. Butter- 

The fading samples genera f a g for a vehicle velocity of 

worth filter are shown in Fig A-7 “*/;f;/„V“uLy is 1-5 GHz and 

20 MPH and 55 that the fading samples for the 

the elevation angle is 90 . It cases, as 

lower speed vary much slower f-eauency wiggle along the slower varymg 
it can be observed, there is a g - out-of-band attenua- 

fading trajectory. This is due enough, so that the 

tion of the third-order Butterworth filt«« 

high-frequency noise p^^ducf some extra phase jitter 

the wiggle. The lugh-frequeu y ^ the simulated system performance 
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To validate the fading simulator, the statistical frequency of a set of 
fading samples is plotted against the theoretical density function. Fig. A. 9 
shows one of these plots for the envelope of the samples in a Rician model 
with K = 10. Similarly, Fig. A.IO shows the phjtse of the fading samples. 
It can be seen that in both cases the histogram agrees very weU with the 
theoreticed probability density function. 

A. 1.2 Simulation of the receiver 

In this section we describe in some detail the general structure of the receiver 
simulation subpackage. 

In particul 2 ir; 

• The received signal may be a binary, quaternary, or octonary CPM 
signal. 

• The front-end filter may have a Gaussian, LSRC, or RC response. 

• The differential detector performs the function ®(t)*a;*(t-T,), where as 
usual T, denotes the symbol interval and * denotes complex conjugate. . 

• The post-detection filter may be Gaussian, LSRC, or RC. 

• The “ph 2 ise measurement” block performs the following operation: 


R[y{t + A)]' 

where, for Q = 2M — 1, 


A - 1 

for full response, 

l» 

for partial response, 

while for Q = M we have 


A - 1 

for fiiU response, 

~ 1 T./2 

for partisJ response. 


A. 2 Eye patterns and their interpretation 

In this section we first briefly describe the significance and the interpretation 
of the eye patterns of the signals obtained at the output of the differenti 2 il 
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detector. Theo «o thaU a nt^ber “g 

“ ‘s^^ofprlfding must be 

CPM, but also for full-response sc eme ^ pulses with truncation length 

contained inaT,-secondm^^^^^^ 

rdt^‘ir:t;et:ltL resold bp verifying their contiatenc, »ith the 
nUTeme t—e “ pC of^a CPM signal he writteo in the form 
0(f,a) = 2irh £ a„q(t-nT,), 
fn 4-1 ±(M - 1)}, at the output of the d^erenti^ 

lel^totfor *»' 

0(t - z, a) : ^ 

^ f* r,T\ 2t/i T a„(i(t - (n + 

^0(t, a) = 2irh an«(t - nZ) - 2irh 

!;‘dm^ifl^T“‘aernd^%S:l^ 

interval (0,iVT,). i,ee v - 1 and consider a partial-response 

let u. take, for iUuatratton’s .^e, N 1, ^d P ^ 

CPM with pvdse length I symbol intervals. Eq. (A.r) can 
rewriting it in the form 

(A.3) 


Ad(t,.) = 2vl. t K-a..i),(<-»T,) + »'mt-b. 

n=-k-L+l 


itT, < 1 < (k + i)r.- 

From the above egnation, -a can obsenie that the eye diagram should 
reflect the shape of the q{-) function. 
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Example 1. 

For a bineiry full-response CPFSK scheme, i.e., 


g(0 = 


0 

t/2Z 

1/2 


t < 0 

0<t<T, , 
t > T, 


we have 

f — hT 

A9{t,a) = wh{ak - ak^i) — — — -+xhak-i, kTs < t < {k + 1)T,. (A.5) 
The corresponding eye diagram is shown in Fig. A.ll. □ 


Example 2. 

Consider now a quatemairy CPM scheme with h = 0.25, and a GMSK pulse 
with truncation length Lx = 1. We have for this case 

= 2irk{ak — ait-i) 9 GMSK(* - kTs) -f irhok-i, kZ <t<{k+ l)Tj. 

(A.6) 

There axe 4 possible initial states, i.e., {— 3 t , —irh,Trh, Sxh}, and four transi- 
tions stemming from each state. For any given state, the slope of the transi- 
tion, i.e., the coefficient of the g(-) function, will depend on the present state 
ak-\Trh and on the input s3rmbol Ofe. For example, the transitions stemming 
from state —Zirh are 

0, iTfhq{t - kZ), ; 6vhq{t - kZ), ; Sirhq(t - kZ)- 
The corresponding eye diagram is shown in Fig. A.12. □ 


Example 3. 

Consider a partial-response binary CPM scheme with rectangular pulses, 
and L = 2. We have in pjirticular 

f 0 t<0 

tl2LZ 0<t<LZ , 

[ 1/2 t> LZ 


g(0 = { 
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Figure A.ll: Eye diagr2im of binary CPFSK. 





Figure A.12: Eye diagram of quaternary CPM with GMSK pulse, h = 0.25, 
Lx = 1 . 


VJ 
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t-kT 

A6{t,a) = 2Trh{ak - 

t ~ kT 

+2irh{ak-i - qfc-2) ' 
+x/iajt_2* 


The corresponding eye diagram is shown in Fig. A.13. □ 
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Appendix B 

Hardware Implementation 


B.l Introduction 

This appendix describes the work done on the firmware implementation of 
the 4800 bps treUis coded CPM modem. The parameters associated with this 
hardware have been determined by the appropriate simulations and analysis. 
The overall functional block diagrams, for the full response and the partial 
response, are shown in Fig. B.l and Fig. B.2 where the transmitter includes 

the following modules: 

• Multiplexer. 

• Trellis encoder. 

• Interleaver. 

• Precoder(not used in the fuU response due to performance degrada- 
tion.) 

• CPM modulator. 

• D/A converter. 

• Lowpass filter. 

• IF frequency converter. 


Accordingly, the receiver includes; 
• Bandpass filters. 
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• Clock recovery. 

• DifFerenti 2 il detector. 

• Deinterleaver. 

• Trellis decoder. 

The end-to-end system requires an extensive digital signal processing, math- 
ematical calculations, and logical operations. The TMS320C25 DSP micro- 
processor (TMS) performs a major part of the mathematical operations. 
The supplemental digital and analog circuits will serve as peripheral hard- 
ware for frequency conversion, filtering, and interfacing. 


B.2 Transmitter 

The transmitter hardware interface dis^am is shown in Fig. B.3. The TMS 
performs, at baseband, the trellis encoding. Interleaving, and CPM modu- 
lation. Additionally, the t imin g for all other circuits m the transmitter is 
coded on TMS. The resulting baseband signals I(t) and Q(t) signals are then 
upconverted to the desired IF level The IF frequency is chosen to be 28.8 
KHz. The transmitter circuits edso needs synchronized timing signals that 
will include: 

• A 4.8 KHz pulse for data source. 

• A 2.4 KHz pulse for symbol rate. 

• A 19.2 KHz pulse for sampling. 

B.2.1 Baseband Processing 

The data, at the rate of 4800 bps, may be provided by a codec or any 
terminal processor. The packetized data is fed to TMS in the form of 2 
parallel bits each at the rate 2.4 Kbps. The data is encoded, in TMS, 
with a rate 2/3 trellis encoder, producing an aggregate rate of 7200 bps. 
Equivalently, the encoder output is generated at the rate of 2400 symbols 
per second (3 bits per symbol). The interleaver block is either (16 by 8) 
or (16 by 16). Implementation of both the trellis encoder and interleaver 
is done by TMS. The interleaved output symbols are input to the CPM 
modulator. The CPM signals are represented by 8 samples per symbol. 
Therefore the required sample rate is 19.2 KHz. This would mean that the 
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TMS must complete the entire entire operations, for each 

1/19200 = 52 Micro seconds. A considerable speed is achieved by the use 

“X TMstft.^eS 

in Fig. B.4. 

B.2.2 IF Converter 

The block dianam of the IF upconvetter it shown in Fig. B.5. Here lk« 

”of287XThe"^^^^^^ 

fiist is a fourth order linear phase lowpass filter for the an og an 
Wdwidth is set at 2.5 KHs. The «cond, on the other 
order chebychev. The transmitter circmt schematics ate shown g. • > 
Fig. B.7, Fig. B.8, and Fig. B.9. 

B.3 Receiver 

The receiver is considerably more compUcated than the ttansnntt«. A M 
“IS receL will he implemented in accortUnc. with the haste block ^ 
gram shown in the Fig. B.IO. Du. to the snbstantial f 

Svolved in the transmitter-receiver operations, a Microseconds 

to accommodate the maximum allowed processing time of 52 Microsecond 
There are eight sets of I and Q signal samples per symbol. The 
^ScSSL; have thus been spHt between four TMS microprocessor systems. 

The assignments are; 

Systems one and two ; perform the upconversions, downconversions , and 

the associated lowpass filtering operations. 

System three ; will do all the baseband processings. 

System four: wiU perform the trellis decoding and deinterleavmg. 
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B.3.1 IF Sampling 

The 28.8 KHz IF signed is analog filtered for noise reduction and then ap- 
plied to the input of the A/D converter. The filtered waveform is sampled 
at 115.2 KHz. The edignment of samples is achieved through a first order 
interpolation of in-phase samples. Here there aue 8 groups of I amd Q sam- 
ples. The groups of the I, Q signads are separated by 52 Microseconds. The 
A/D output I auid Q samples aue separated by 8.681 Microseconds. This 
separation sets the upperbound on the converter speed. Therefore, a fairly 
high speed converter is required to meet the required tumng . To retain 
the baseband component lowpaiss filtering is carried out by using TMS im- 
plemented appropriate FIR filters. Fig. B.IO illustrates the basic receiver 
architecture at the IF level. 

The ground noise is expected to be larger than 1/256 of the input signal. 
Therefore, an 8 bit resolution associated with the A/D device should provide 
sufficient accuracy. The schematics of TMS systems one and two are shown 
in Fig. B.ll and Fig. B.12. 

B.3.2 Baiseband Operations 

The baisebamd operations include filtering, clock recovery, amd differential 
detection. These functions are implemented on the third TMS system. The 
block diagram of the differentiad detector is shown in Fig. B.13. Fig. B.14 
amd Fig. B.15 illustrate detauls of the third amd fourth TMS systems. 

B.3.3 Decoder and Deinterleaver 

The fourth TMS system will perform operations of deinterleaving amd Viterbi 
decoding. The schematics of this system is shown in The TMS is interrupted 
at the rate of 2.4 KHz. After the filtering, Doppler correction, amd demodu- 
lation the data will be deinterleaved amd decoded. Then the decoded symbols 
are loaded into the shift registers. 


B.4 Status of the Hardware 

1. The tramsmitter hardware has been completed and tested. The circuits 
include : 

• -The TMS microprocessor system. 

• Peripherad circuits. 
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• IF up conversion circuits. 

The TMS softwaxe for the transmitter is complete. It does not include 
the encoder. The eye diagrams corresponding the CPM signal with 
IRC pulse, have been photographed and are shown in Fig. B.16 and 
Fig. B.17. The current system runs on a 2400 bps data rate. Thus 
transmitter hardware must be upgraded to accommodate the required 
4800 bps data rate. This transmitter box, with operating instructions, 
has been delivered along with the final report. 

2. As for the receiver, a total of four TMS microprocessor systems make 
up the receiver. The hardware for these systems have been completed. 
The components include: 

• The memory system. 

• Memory peripheral circuits. 

• System’s timing, 

• A/D converter. 

• Interrupt circuits. 

• Handshaking logics for interfacing of the multiprocessor opera- 
tions. 

Moreover, except the decoder, the software for the baseband receiver 
has been completed. In summary, additional time is required to com- 
plete the hardware implementation of the differential CPM system. 
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The Partial response system 


Figure B.2 


Channel 
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Transmitler block diagram 


Figure B.3 
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Figure B.4 
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Figure B.5 
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Figure B.6 
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Figure B.7 
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Figuie B.8 
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Figure B.9 
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Receiver Architecture at IF band 


Figure B,10 






appendix b. hardware implementation 


274 



Board one : A/D, LPF, Detector ( I component) 


Figure B.ll 
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Figure B.12 
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The full response system 


Figure B.13 
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Figiire B.14 
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Board 4 ; decoder and deinterleaver 


Figure B.15 
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Figure B.16 
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